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Abstract: Gastropods have a dense fossil record dating back to the Late Cambrian. Intuitively, this would appear to aid phylogenetic reconstructions. 
However, workers question both whether gastropod shell characters are phylogenetically informative and whether stratigraphic data can be used to test phy¬ 
logenetic hypotheses. Both questions are addressed with an analysis of 82 species of the Lophospiroidea (= Lophospiridae + Trochonematidae) from the 
Ordovician and Silurian. Compatibility analyses of 95 shell characters shows that characters are far more compatible than one would expect from homo- 
plasy-saturated data. However, compatibility among the 73 characters that vary among the earliest lophospirids decreases over time, which suggests that 
later species introduced homoplasy. Simulations using 95 characters and sampling similar to that observed for lophospirids show that parsimony performs 
poorly relative to methods that incorporate stratigraphic data such as stratocladistics. An alternative approach is used here. The first step is to estimate the 
likelihood of a hypothesized tree given observed character congruence (/. e. parsimony length) using simulations. The second step uses two different statisti¬ 
cal tests to estimate the likelihood of hypothesized trees given observed stratigraphic data. Likelihoods then are combined to evaluate trees. The resulting 
likelihood tree is nearly 30 steps longer (378.4 versus 350.5), but is considered more likely (given a 350.5-step matrix) than a tree of 350.5 steps. Both trees 
suggest that budding cladogenesis (where ancestors co-exist with descendants) was the most common pattern of speciation, although the likelihood tree is 
more emphatic on this point. Both suggest a trend towards increasing numbers of ornate species; however, whereas parsimony suggests the differential 
diversification of an ornate clade, likelihood suggests a strong tendency for inornate ancestors to have ornate descendants. A genus-level taxonomic revision 
is provided that is consistent with both trees and that attempts to reflect historical diversity patterns. 
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Gastropods have a rich fossil record coupled with 
extensive ecologic and morphologic diversity among extant 
species. These factors make the clade an excellent model 
for testing macroevolutionary hypotheses (Bieler, 1992). 
Most macroevolutionary hypotheses make predictions 
about phylogenetic patterns. Therefore, robust phylogenet¬ 
ic estimates offer means of testing alternative macroevolu¬ 
tionary hypotheses. Examples of such hypotheses that have 
been tested in a phylogenetic context with fossil gastropod 
data include hypotheses about speciation patterns (Wagner 
and Erwin, 1995), predator-prey escalation (Carlson and 
Vermeij, 1996), phylogenetic constraints (Wagner, 1995a), 
and long-term trends (Wagner, 1996). 

Unfortunately, fossilized character data for gas¬ 
tropods are limited almost entirely to shell characters, 
which likely are highly homoplastic ( e . g . Harasewych, 
1984; Kool, 1993). Phylogenetic methods such as parsimo¬ 
ny likely will perform poorly with shell data because much 
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character congruence will reflect convergence rather than 
shared ancestry (Felsenstein, 1984; Archie, 1996). The 
dense fossil record of gastropods offers a potential antidote 
to the homoplasy problem because phylogenetic hypotheses 
make necessary predictions about durations, which in turn 
make probabilistic predictions about stratigraphic ranges. 
Accordingly, several workers have used stratigraphic data 
to test whether characters linking taxa are better explained 
as convergences than as homologies ( e . g. Fisher, 1991, 
1994; Huelsenbeck, 1994; Cheetham and Jackson, 1995; 
Wagner, 1995b). 

This paper will attempt to address whether shell 
character data retain phylogenetic signal and whether incor¬ 
porating stratigraphic data improves estimates of phyloge- 
nies. Shell character and stratigraphic data for Ordovician 
and Silurian members of the Lophospiroidea (= 
Lophospiridae + Trochonematidae) are analyzed using 
compatibility, parsimony, and maximum likelihood tech¬ 
niques to contrast alternative hypotheses about lophospiroid 
relationships and character evolution. The implications of 
the different results for macroevolutionary hypotheses and 
lophospiroid classification also are discussed. 
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DATA AND ANALYSES 

SHELL CHARACTER DATA, ANALYZED SPECIES, 
AND PARSIMONY 

A previous phylogenetic analysis of lophospiroids 
used 57 characters and 150 character states to estimate rela¬ 
tionships among 55 Ordovician species (Wagner, 1995b). 
The analysis presented here uses 95 characters encompass¬ 
ing 257 character states for 82 Ordovician and Silurian 
species. The increased number of character states reflects 
in part the addition of Silurian species with characters and 
character states unobserved among Ordovician species. 
Also, the characters were re-coded following Wagner [in 
press (a)]. Many characters from the 1995 analysis are 
divided into two or more characters that better reflect varia¬ 
tions observed among Early Paleozoic gastropods. 
Nevertheless, 95 characters might appear to be a large num¬ 
ber for a group of gastropods that previously were divided 
into only two families. However, a survey of numerous 
phylogenetic studies of gastropods shows that the ratio of 
shell characters to analyzed taxa increases as the taxonomic 
level of the analysis becomes finer [Wagner, in press (a)]. 
Also, lophospiroids are a morphologically diverse clade. 
Many more characters are needed to described this diversity 
than are needed to describe any one species. 

The characters are listed and described in Appendix 
1. The character matrix is given in Appendix 2. 
Continuous characters are divided into ordered series using 
segment coding and then de-weighted so that the maximum 
difference equaled one step (Chappill, 1989). As a result, 
alternative trees often have fractional lengths. Other multi¬ 
state characters are treated as ordered if they represented a 
logical geometric series; otherwise, they are considered 
unordered. Some characters, such as the shapes and dimen¬ 
sions of the right and left ramps, vary both in symmetry 
(/. e. left and right being nearly identical) and dimensions. 
For such characters, symmetry versus asymmetry is coded 
as one character. Left and right characters are coded sepa¬ 
rately if they varied independently among asymmetric 
species. However, left and right characters necessarily 
covary on symmetric species, which leads to a conundrum: 
devising a coding scheme where the difference in left and 
right ramp shape accounts for one difference between two 
symmetrical species, two differences between two asym¬ 
metrical species, and two or three differences between a 
symmetrical and an asymmetrical species (one difference 
reflecting symmetry.asymmetry, with the third difference 
apparent only if both the left and right ramps differ). Step 
matrices (Swofford and Olsen, 1990) offer a potential solu¬ 
tion, but these result in exorbitant run times. Instead, asym¬ 
metrical characters are weighted one-half of 
presence/absence characters, making changes among sym¬ 
metrical species equal to one step. Trees are rescaled using 


a separate program that insures no “chimera” reconstruc¬ 
tions ( e . g. hypothesized symmetrical ancestors with differ¬ 
ent left and right sides), with the tree lengths now scaled as 
they would be by a step matrix. 

Broader phylogenetic analyses [Wagner, in press 
(a)] suggest that lophospiroids are nested within the 
Murchisoniina, with species assigned to the genus 
Ectomaria representing the immediate outgroup. 
Therefore, Ectomaria adelina, which is both the earliest 
known Ectomaria and contemporaneous with the oldest 
lophospiroids, is included as an outgroup. Another contem¬ 
poraneous murchisoniinae species, Hormotoma simulatrix, 
also is included as an outgroup. 

Heuristic parsimony analyses using PAUP 4.0* 
(Swofford, 1998) found 36 trees of 350.8 steps. The illus¬ 
trated tree (Fig. 1) is the parsimony tree with the fewest 
stratigraphic gaps, which satisfies the secondary role for 
stratigraphy advocated by Smith (1994). [This is equiva¬ 
lent to Clyde and Fisher’s (1997) “Analysis 2.”] The tree is 
discussed in greater detail below, but it is presented here to 
introduce the potential errors for the following sections. 
The 50 random addition replicates found an additional six 
islands of trees (Maddison, 1991) with minimum length 
trees of 380.5 or fewer steps. The phylogenetic reconstruc¬ 
tion (Fig. 1) accounts for both stratigraphic ranges and 
character distributions. Both are important because species 
identical to reconstructed ancestral morphologies obviate 
the need for hypothesized range extensions (Fisher, 1991; 
Smith, 1994). 

DO SHELL CHARACTERS CONTAIN 
PHYLOGENETIC SIGNAL? 

Bandel and Geldmacher (1996) criticized a previous 
phylogenetic analysis of Ordovician lophospiroids by 
Wagner (1995b) because the resultant phylogenetic esti¬ 
mate linked Trochonema and allies to lophospirids. Bandel 
and Geldmacher (1996) considered these to be distant rela¬ 
tives based on Triassic specimens. Those authors attributed 
Wagner’s (1995b) results to pervasive homoplasy among 
shell characters. There is some support for Bandel and 
Geldmacher’s (1996; hereafter B&G) hypothesis. 
Parsimony reconstructions based on soft anatomy suggest 
high frequencies of homoplasy among shell characters 
( e . g. Carlson and Vermeij, 1996; Haasl, 1997). The utility 
of shell characters (on which this analysis necessarily 
relies) must be questioned and the B&G hypothesis must be 
tested before proceeding with further analyses. 

Hypothesized phylogenies predict hierarchical 
structure among homologies. Conversely, a hypothesis of 
pervasive homoplasy (such as the B&G hypothesis) pre¬ 
dicts little hierarchical signal among the characters. 
However, the B&G hypothesis does not predict an absence 
of congruence. Although the terms often are confounded 
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Fig. 1. Phylogenetic interpretation of one of 36 parsimony cladograms. This tree implies less stratigraphic debt than does any of the other 35 trees. Thin 
solid lines show estimated phylogenetic connections. Dashed lines show hypothesized unsampled lineages and taxa that are implicit to the phylogenetic 
hypothesis (j. e. range extensions sensu Smith, 1988, 1994; = ghost lineages and taxa sensu Norell, 1993). The tree posits 350.5 steps [consistency index (Cl) 
= 0.316; retention index (RI) = 0.807) and 102 units of stratigraphic debt, with debt units based on the substages marked on the time scale. Species identical 
to hypothesized ancestors are linked directly to nodes and are considered ancestral to their apomorphic sister taxa here. [Trochonemella “ G . pulchellum” is 
identical to Trochonemella BMNH 36364 in Appendix 2.] 


(e. g. Bryant, 1992; Carpenter, 1992), hierarchy and congru¬ 
ence are not synonyms: congruence exists in non-hierarchi- 
cal matrices (Kallersjo et al , 1992; Alroy, 1994) and even 
in randomly generated matrices (Archie, 1989; Faith and 
Cranston, 1991). However, whereas hypotheses of phyloge- 
ny and homology predict hierarchical congruence among 
characters, hypotheses of homoplasy predict non-hierarchi- 
cal congruence. (The one exception to this expectation is 
when multicharacter complexes are homoplastic and pat¬ 
terned so that character combinations appear in the same 
order multiple times; however, this begins to violate another 


important assumption of phylogenetic methods, i. e. the 
independence of characters.) 

Several workers ( e . g. Meacham, 1984, 1994; 
Sharkey, 1989, 1994; Alroy, 1994) have suggested using 
character compatibility (Le Quesne, 1969; Estabrook et al ., 
1975) to evaluate hierarchical content of matrices. Two 
binary characters with pairings of (00), (10), (01), and (11) 
must have homoplasy in one or both characters. Such 
incompatible characters necessarily lack global hierarchy. 
Two binary characters with only three of the possible pairs 
[e. g. (00), (10), and (01) or (11)] are considered compatible 
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because they can be reconstructed on some trees without 
homoplasy. Characters with little homoplasy should form 
many compatible pairs whereas those with numerous 
homoplasies should form many incompatible pairs. 

Alroy’s (1994) Permutation Compatibility (PC) test 
is operationally similar to the Permutation Tail Probability 
test (Faith, 1991). PC tests whether the hierarchical struc¬ 
ture of a character matrix is significantly greater than that 
of a randomly produced matrix. Alroy (1994) divided com¬ 
patibility into two types: hierarchical [= “direct” of Sharkey 
(1994)] and general (“indirect”). The former describes 
pairs where combinations show a hierarchical arrangement 
[e . g. (00), (01), and (11)]. The latter describes pairs where 
there is neither necessary homoplasy nor implied hierarchy 
[e. g. (00), (01), and (10)]. Because phylogenetic hypothe¬ 
ses explicitly predict hierarchy, I used the PC test with 
1,000 randomized matrices to determine if randomly per¬ 
muted matrices ever retained the same level of hierarchical 
compatibility as did the real matrix. Lophospiroid character 
data yielded 1,389 hierarchically compatible pairs out of 
8,930 possible pairs. (The test has been amended here to 
compare characters, not character states, as was done in 
Alroy’s implementation). Permuted matrices typically 
showed approximately half as many pairs (Fig. 2). Thus, 
the matrix has more hierarchical structure than expected if 
congruence largely reflected random homoplasy, allowing 
us to reject the B&G hypothesis. 



Hierarchically Compatible 
Character Pairs 

Fig. 2. Results of the Permutation Compatibility test (Alroy, 1994). The 
arrow denotes the observed value whereas the histogram gives the values 
generated by 1,000 random permutations of the character data. 
Hypothesized phylogeny predicts hierarchical compatibility whereas other 
hypotheses of character congruence (e. g. general convergence) do not. 
Thus, these results support (but do not demonstrate) the idea that there is 
phylogenetic signal among these characters. 


DO WE NEED STRATIGRAPHIC DATA TO 
IMPROVE OUR ESTIMATES OF LOPHOSPIROID 
PHYLOGENY? 

The presence of some hierarchical structure is no 
guarantee that there is not sufficient homoplasy to mislead 
parsimony. Simulation studies show that even low rates of 
randomly accrued homoplasy will result in inaccurate par¬ 
simony reconstructions of phylogeny ( e . g. Mooers et al ., 
1995). Finite numbers of character states (see Wagner, 
1998c) insure that the probability of chance homoplasies 
confounding parsimony is not infinitesimal (Felsenstein, 
1978). Patterned homoplasies ( e . g. functional complexes) 
violate assumptions of character independence and also 
lead to inaccurate parsimony reconstructions (Lamboy, 
1994; Archie, 1996). Finally, rates of character change that 
are not adequately reflected by character weighting (see 
Felsenstein, 1981) also yield inaccurate parsimony trees 
(Kuhner and Felsenstein, 1994). 

Neontologists can test phylogenetic hypotheses 
derived from one data set by examining different character 
sets ( e . g. different types of molecules). The hypothesis 
that congruence from the original data set reflects phyloge¬ 
ny predicts very similar patterns of congruence in other 
data sets (assuming that changes in one character set do not 
affect changes in the other; Templeton, 1983). 
Unfortunately, only a single character set is available for 
extinct taxa such as lophospiroids. However, phylogenetic 
hypotheses also make necessary predictions about temporal 
durations of taxa (see, e. g. Smith, 1988). Predictions about 
durations make probabilistic predictions about stratigraphic 
ranges (Paul, 1982; Strauss and Sadler, 1989). If a hypoth¬ 
esis makes predictions (necessary or probabilistic) about a 
data set, then those data offer a test of the hypothesis. 

Additional factors suggest that stratigraphic data 
might improve estimates of lophospiroid phylogeny. 
Simulation studies [Wagner, in press (b)] show that phylo¬ 
genetic error exaggerates the phylogenetically implied 
range extensions (Smith, 1988, 1994; “ghost” lineages and 
taxa sensu Norell, 1993) far more often than error underes¬ 
timates range extensions. The exaggeration becomes worse 
as frequencies of change per character per branch increase, 
but it is pronounced at the frequencies of change posited by 
the parsimony tree for lophospiroids (/. e. f = 0.042 per 
character per branch). Thus, the long gaps in sampling 
posited in Fig. 1 are evidence that synapomorphies (sensu 
Sober, 1988) are homoplastic rather than homologous. For 
lophospirids, most gaps are highly improbable given 
observed sampling distributions (Wagner, 1995b). 

Another concern is that the hierarchical signal of the 
character matrix might decrease as geologically younger 
taxa are added. Suppose that geologically older species 
have morphotypes (00), (01), and (01), but that later species 
introduce morphotype (10). The characters no longer are 
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Number of Species 

Fig. 3. Effects of geologically younger taxa on overall compatibility. As 
taxa are added over time, compatibility decreases. The white line gives the 
compatibility among all characters after X species are added to the matrix. 
Crosses give the compatibility among only the 73 characters informative 
for “early” (/. e. pre-Caradoc) relations. Thin black line gives the compati¬ 
bilities for all characters when stratigraphic ranges are assigned at random. 

compatible, and the younger species necessarily imply 
reversal or convergence. Conversely, suppose that the older 
species had states (00) and (11) and some geologically- 
younger species had (10). This suggests that the “younger” 
species is descended from a phylogenetic intermediate that 
was not sampled initially. The former scenario predicts 
decreased character compatibility whereas the latter does 
not. 

I examined the effect of geologically younger taxa 
on character compatibility among lophospirids using all 
character and also only the 73 characters that vary among 
early lophospiroid species (Fig. 3). (Late-appearing char¬ 
acters are not relevant to relationships among the oldest 
taxa.) Compatibility decreases markedly over time. 
Sampling so poor that species were effectively sampled at 
random through time relative to their phylogenetic position 
(as implied by abundant ghost lineages) predicts a very dif¬ 
ferent pattern. One thousand randomizations of stratigraph¬ 
ic ranges show that compatibility should begin lower than 
observed and that the decrease over time should be much 
less precipitous (Fig. 3). The observed patterns are consis¬ 
tent with the idea that Late Ordovician and Silurian lophos- 
piroids exhibited homoplasy among the characters that dis¬ 
tinguish Early-Middle Ordovician lophospiroids. 

A final concern is whether methods incorporating 
stratigraphic data perform better than does parsimony. 
(“Parsimony” here and throughout this paper denotes mini¬ 
mum steps evolution [Edwards and Cavalli-Sforza, 1964; 


Kluge and Farris, 1969]; other methods using parsimony 
criteria \e. g. stratocladistics] are labeled differently.) To 
test this, I simulated morphologic evolution using the num¬ 
ber of characters and states per character apparently avail¬ 
able to lophospiroids. Character evolution was ordered or 
unordered reflecting the assumptions of the phylogenetic 
analysis (Appendix 1). In one round of analyses, log-prob¬ 
abilities of character change were proportional to the rela¬ 
tive weighting of characters (i. e. “equiprobable change”; 
see Felsenstein, 1981). In a second round, rates of change 
varied among characters at random (/. e. “variable 
change”). Two different speciation models were used, one 
in which ancestral morphotypes survived cladogenesis and 
could produce any number of descendants (“budding clado¬ 
genesis”) and one in which ancestral morphotypes became 
pseudo-extinct while giving rise to two distinct descendants 
(“bifurcating cladogenesis”). Sampling parameters derived 
from lophospiroid data (see Appendix 3) were used to sam¬ 
ple six species. (A six-species limit was imposed so that 
exhaustive searches could be used.) Parsimony analyses 
then were run using the simulated character matrices. 
Matrices were maintained only if parsimony estimated the 
same steps per sampled taxon as the real lophospiroid tree 
(/. e. 350.5 steps per 82 taxa ? 25.5 steps per six taxa; here¬ 
after: a 25.5-step matrix). 

The simulations then had character matrices and 
fossil records that were similar to those observed for 
lophospiroids. The stratigraphic data then were used for 
three additional phylogenetic analyses: reweighting with 
the Stratigraphic Consistency Index (SCI; Huelsenbeck, 
1994), stratocladistics (Fisher, 1994), and sieving with 95% 
confidence intervals (CIS; Wagner, 1995b). The results of 
the phylogenetic analyses were contrasted using Robinson 
and Foulds’ (1981) metric, which measures the proportion 
of nodes that agree on two trees. (If multiple optimal trees 
were found, then the mean error was used.) 

All methods incorporating stratigraphic data outper¬ 
formed parsimony when simulated character matrices and 
simulated sampling matched observed parameters for 
lophospiroids (Fig. 4). In the best case, stratocladistics cor¬ 
rectly reconstructed over 50% of the cladograms. 
Parsimony, meanwhile, correctly reconstructed under 15% 
of the cladograms. Stratigraphic methods are less success¬ 
ful when rates of character evolution are not proportional to 
character weighting, something that is true of parsimony 
and other phylogenetic methods (Kuhner and Felsenstein, 
1994). However, stratigraphic methods are slightly less 
affected by such variation than is parsimony. 

In summary, lophospiroid character data show too 
much compatibility to be dismissed as noise. However, 
reasons to question the parsimony estimate of phylogeny 
include: (1) many implied but statistically improbable sam¬ 
pling gaps, (2) decreasing hierarchical signal over time, and 
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Fig. 4. Success of four phylogenetic methods for six-taxon simulations. 
The numbers reflect the number of incorrect nodes (“0” indicates a correct 
tree). The simulations use 95 characters and the same number of possible 
character states and character state orderings per character as lophospiroid 
data. Sampling intensities mimic observed sampling. “Equiprobable” evo¬ 
lution used probabilities of character change proportional to the character 
weight, whereas “variable” evolution used randomly assigned probabilities 
of change. Results are from budding cladogenesis, in which ancestral 
species can survive to yield any number of descendants, but the simula¬ 
tions produce nearly identical results when using bifurcating cladogenesis. 
Note that all three methods using stratigraphic data outperform parsimony. 
See text for additional discussion. 


(3) superior performance of methods incorporating strati¬ 
graphic data in simulations. 

WHY MAXIMUM LIKELIHOOD APPROACHES 
ARE USEFUL 

How best to incorporate stratigraphic data into phy¬ 
logenetic analyses is a non-trivial issue. I reject the posi¬ 
tion that stratigraphic data cannot improve phylogenetic 
estimates because such data are not hierarchically distrib¬ 
uted ( e . g. Eldredge and Cracraft, 1980; Rieppel and 
Grande, 1994; Smith, 1994). This only means that strati¬ 
graphic data do not offer inductive statements about phy- 


logeny. If the goal is to test the predictions of phylogenetic 
statements, then one criterion applies: do estimated phylo- 
genies make predictions about stratigraphic data? The syl¬ 
logism is straightforward: if an inferred phylogeny is accu¬ 
rate, then taxa originated by particular times; if taxa origi¬ 
nated by a particular time, then there is a some probability 
that the known fossil record would be observed. Of course, 
the “deduction” is necessarily fuzzy, as stratigraphic data 
cannot demonstrate that a lineage had originated by any 
point in time, only that it is highly unlikely to have done so 
given observed data. 

Unfortunately, the methods incorporating strati¬ 
graphic data used above all are unsatisfactory. Reweighting 
trees with SCI suffers because the SCI is greatest on pecti¬ 
nate trees (Siddall, 1996). Sieving trees with confidence 
intervals is not obviously predisposed toward favoring par¬ 
ticular topologies. However, CIS suffers from focusing on 
particular nodes rather than the whole tree (Wagner, 
1998b). Multiple independent tests and a traditional signif¬ 
icance value of 0.05 should yield Type I errors (/. e. incor¬ 
rectly rejecting the null) one time in every 20 tests. More 
damning, multiple tests should yield Type II errors (i. e. 
incorrectly accepting the null) even more frequently. For 
example, if there are ten gaps that are significant at a p- 
value of 0.10 each, then the CIS method will accept all of 
them. However, we expect only one of those gaps to be 
real. 

An alternative approach with sounder logical justifi¬ 
cations is to use likelihood to evaluate hypothesized trees. 
Stratocladistics actually represents a simple likelihood test. 
Summing steps is equivalent to summing negative logs of 
character change probabilities (Felsenstein, 1981) if: (1) 
probabilities for each character are the same among all 
clade members, (2) characters never change twice on the 
same branch, and (3) characters always evolve indepen¬ 
dently. The most likely tree then is the one invoking the 
fewest steps. Summing the negative log-probabilities of 
implied gaps (/. e. stratigraphic debt setisu Fisher, 1991, 
1994) estimates the log-probability of a tree. Again, the 
most likely tree is the one implying the fewest stratigraphic 
gaps. The log-likelihood of the entire tree now equals the 
sum of stratigraphic debt and morphologic steps if the 
stratigraphic debt is weighted as: 

-In (1 - R) = W *-In (P [c]) 

where 1 - R is the probability of not sampling a taxon over 
a particular interval (Foote and Raup, 1996), W is the 
weight of stratigraphic debt relative to characters, and P [c] 
is the probability of character change along a branch. 

Stratocladistics will retain nodes implying improba¬ 
ble gaps if they are supported by numerous synapomor- 
phies. This offers a control for Type I errors which is lack¬ 
ing in the confidence interval method. Clyde and Fisher 
(1997) further justified stratocladistics on the grounds that 
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it uses the same optimality criterion to evaluate all of the 
data, which is the justification of “total evidence” analyses 
( e . g. Kluge and Wolf, 1993). However, stratocladistics 
requires a priori assumptions about sampling intensity ( R) 
and frequencies of character change. As unknowns, both 
should be tested rather than assumed. Fortunately, the 
stratigraphic debt for any one tree implies an /?, the likeli¬ 
hood of which can be assessed analytically (Foote, 1997). 
R's predicting some debt are more likely than R’s predict¬ 
ing no debt (Wagner, 1998a), which means that statistically 
optimal trees can be “suboptimal” by stratocladistic criteria. 

Nevertheless, the basic format of stratocladistics 
serves as a useful logical template. Because stratigraphy and 
morphology are independent data sets, the overall likelihood 
of any hypothesized tree is the likelihood of the tree given 
the character data times the likelihood of the tree given the 
stratigraphic data (Edwards, 1992). This is the approach that 


will be taken in this paper (see also Wagner, 1998a). 

CALCULATING LIKELIHOOD OF A HYPOTHE¬ 
SIZED TREE GIVEN THE CHARACTER DATA 

Congruence is an observable datum of a character 
matrix that is summarized by the parsimony length. The 
trees of interest are treated as hypotheses to be tested here. 
The likelihood of a hypothesized tree length (HL) yielding 
a parsimony length (PL) is proportional to the probability 
of observing the datum (PL) given the hypothesis (true tree 
length) (Edwards, 1992). Thus: 

L [HL I PL] ocP [PL I HL] 

P [PL I HL] cannot be calculated analytically, but 
simulations can be used to determine the frequencies of 
PL’s derived from trees of known lengths. The simulations 
illustrated in Fig. 4 were based on six-taxon trees with 
25.5-step character matrices. (Again, six-taxon trees 
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Fig. 5. The distributions of parsimony lengths given known lengths, based on simulations for six taxa with lophospiroid characters. The frequencies of matri¬ 
ces yielding parsimony lengths of 25.5 steps are highlighted. These match the number of steps per taxon as found for all 82 lophospiroids. 
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allowed use of exhaustive searches.) The frequencies of 
PL’s given known lengths (/. e. f [PL ! Known Length] in 
Fig. 5) indicate that a tree of 27 or 27.25 steps has a much 
higher probability of yielding a 25.5-step character matrix 
than does a tree of 25.5 steps. Thus: 

L [HL = 27.25 I PL = 25.5] > L [HL = 25.5 I PL = 25.5] 
This is best appreciated by examining the relative heights 
of the solid gray bars in Fig. 5, which approximate the like¬ 
lihood of HL’s given the observed congruence. The bar at 
27.25 is the highest, indicating that this is the length most 
likely to have yielded the observed congruence. 

One point not discussed above is that the simula¬ 
tions estimate P [PL I KL] only if a particular hypothesis of 
character evolution is assumed. Fig. 5 illustrates the results 
when simulated character evolution is derived from the 
weighting scheme (see Felsenstein, 1981). Deviations from 
this models resulted in longer trees lengths becoming more 
likely and the shortest tree length becoming less likely 
(Wagner, 1998a). Thus, these results represent the highest 
possible likelihoods for parsimony hypotheses. 

We need to evaluate 82-taxon phylogenies to evalu¬ 
ate lophospiroid phylogenies. This is accomplished by ran¬ 
domly selecting 13 six-taxon clades and one four-taxon 
clade (derived from a separate set of simulations not illus¬ 
trated here). Parsimony lengths among random six-taxon 
clusters within lophospiroids varied from 15.8 to 29.6 steps. 
Therefore, true tree lengths were randomly selected from 
the range of true lengths that yielded parsimony lengths of 
15.8 to 29.6 steps (/. e. 15.8 to 54.1 steps). Note that the 
use of a range of true lengths means that estimated proba¬ 
bility does not assume clock-like rates of morphologic 
change (contra Norell and Novacek, 1997). 

Parametric bootstrapping (see Huelsenbeck et cil , 
1996) then assigned parsimony lengths to each true length. 
Assignments used distributions such as those illustrated in 
Fig. 5. The true lengths and parsimony lengths of the 14 
subclades are then summed, which yields the likelihood 
distribution for hypothesized tree lengths given 82 taxa and 
lophospiroid character data (Fig. 6). Log-likelihoods are 
illustrated because the likelihoods of shorter trees are too 
low to be visible on a histogram. 

Parametric bootstrapping never reconstructed any 
trees with both true and parsimony lengths of 350.5 steps. 
(This is not surprising, as the probability of doing so is 
approximately KH 4 , which requires several orders of mag¬ 
nitudes more replications than performed.) Because the 
likelihoods of tree lengths under 356 steps could not be 
estimated from parametric bootstrapping, I assigned log- 
likelihoods of that length (-10.415) to all shorter lengths. 

Alroy (pers. comm., 1998) suggested that the likeli¬ 
hood estimates for the parsimony length might be too low 
because heuristic searches were used to estimate the initial 
parsimony length of 350.5 steps. However, for a hypothe- 



Hypothesized Length (HL) 

Fig. 6. Log likelihoods of hypothesized tree lengths given a 350.5 step 
matrix for 82 taxa with lophospiroid characters. Note that trees 20-30 steps 
longer than the parsimony length are orders of magnitude more likely than 
are trees of 350.5 steps. 

sized length of 350.5 to be as likely as a tree length of 370, 
the actual parsimony length must be around 330 steps. It 
seems improbable that heuristic searches could have missed 
by that much. Other aspects of the test overestimate the 
likelihoods of parsimony lengths (Wagner, 1998a). The test 
does not allow parsimony to confound homoplasy among 
subclades. Also, it uses the best-case model of character 
evolution (see above). All other models result in a lower 
probability of parsimony lengths equaling hypothesized 
lengths. These factors, coupled with the deliberate overes¬ 
timate of parsimony length likelihood described above, 
mean that the tests used herein are quite conservative if one 
wishes to treat the parsimony tree as a null hypothesis. 

CALCULATING LIKELIHOOD OF A HYPOTHE- 
SIZED TREE GIVEN STRATIGRAPHIC DATA 
Continuous Stratigraphic Data 

Wagner (1995b) used continuous stratigraphic data 
to estimate confidence intervals for lophospiroids. 
Huelsenbeck and Rannala’s (1997) maximum likelihood 
test was applied using these data (see Appendix 3). As in 
the 1995 analysis, the “time” scale used here is actually the 
number of sampled units (“horizons”) per stratigraphic 
interval. Thus, gaps through poorly sampled intervals are 
more likely than are gaps through well-sampled intervals, 
even if both intervals are of the same duration. 

Huelsenbeck and Rannala’s (1997) test is modified 
slightly here. The most important modification is that the 
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phylogeny is no longer considered a predictor of species’ 
extinctions. This would be true if species evolved anage- 
netically. However, budding cladogenesis patterns predict 
only the latest possible origins of species. The parsimony 
tree (Fig. 1) illustrates possible examples of budding clado¬ 
genesis, as parsimony finds several plesiomorphic species 
known after their apomorphic sister taxa appear. 
(“Plesiomorphic species” are those that match inferred 
ancestral morphotypes.) Previous phylogenetic analyses 
suggest that budding cladogenesis was the predominant 
speciation pattern among lophospiroids (Wagner and 
Erwin, 1995). The exact probability of data given a 
hypothesized origination now is: 

P [FKA, LKA, H I HFA, /] = f ((LKA - FKA)H - 2 fH e -/ 
(LKA - HFA), (H - 2)! ) when H _ 2 

= /e-/(LKA-HFA) when H = 1 

where FKA is the first-known appearance, LKA is the last- 
known appearance, HFA is the hypothesized first appear¬ 
ance, H is the number of horizons from which the species is 
known, and / is the proportion of horizons within a species 
range from which it is sampled (amended from 
Huelsenbeck and Rannala, 1997: equations 2 and 4). 

The likelihood of a hypothesized origination time 
implicit to the inferred phylogeny now is: 

L [HFA, X I FKA, LKA, H ] = c * IPR(i = U = nodes + 
OTUs-1 , P (FKA, LKA, H I HFA, X) ) 

(amended from Huelsenbeck and Rannala, 1997: equation 5) 
where c is an arbitrary constant that rescales each likelihood 
(L), so that the maximum likelihood is 1.0 (Edwards, 1992). 

Discrete Stratigraphic Data 

The stratigraphic debt implicit to an inferred phy¬ 


logeny implies a particular R (Wagner, 1998a), the likeli¬ 
hood of which can be evaluated (Foote, 1997). (“Sampling 
intensity” here denotes the average proportion of species 
per interval that are sampled; Foote and Raup, 1996). If a 
species invokes four units of stratigraphic debt, then the 
phylogeny implies that the species was sampled in its fifth 
unit after diverging from its closest sampled relative. The 
phylogeny therefore implies a sampling intensity of 0.2 
(/. e. one find in five tries). For an entire phylogeny, 90 
units of stratigraphic debt distributed among ten species 
suggests an average sampling intensity of 0.1. Because 
phylogenies infer latest necessary divergences (Smith, 
1988), stratigraphic debt gives the minimum implied gaps. 
A simple analytic estimate therefore is biased toward over¬ 
estimating R, especially on large phylogenies (Wagner, 
1998a), so R is best estimated using simulations. 

Stratigraphic debt was calculated using the sub¬ 
stages marked in Fig. 1. For the parsimony tree, there are 
102 units of stratigraphic debt. Stratigraphic debts of 102 
for 82 species (two of which are present in the earliest sub¬ 
stage) imply an R = 0.33. (Note that an analytic R = f (80, 
182) = 0.44 is substantially greater.) 

Taxon ranges reflect both R and durations (which 
reflects extinction intensity, fj) (Sepkoski, 1975; Foote and 
Raup, 1996). As sampling becomes worse, ranges will 
decrease regardless of actual durations. In the extreme case 
of nearly infinitesimal R, the few observed species will each 
be known from single localities and all species will have 
ranges of one unit. Such clearly is not the case for lophos¬ 
piroids (Fig. 7A). The most likely hypothesis for observed 
lophospiroid ranges is R = 0.68 per substage (Fig. 7B, with fj 
= 0.52 per substage). Given the R associated with given 




r> 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
D Sampling Intensity (R) 


Fig. 7. Estimates of the likelihood of sampling intensities ( R ) and extinction intensities (/i). A. Distributions of observed ranges, with ranges measured by the 
major divisions shown in Fig. 1. B. Log-likelihoods for hypotheses that particular combinations of fj and R yielded the distribution shown in A. The maxi¬ 
mum likelihood estimate is R = 0.68 per substage and fj = 0.52 per substage. 
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amounts of stratigraphic debt, it now is possible to assess the 
likelihood of any given amount of stratigraphic debt (Fig. 8). 

RESULTS 

BASIC ASPECTS OF THE ALTERNATIVE TREES 

Likelihoods were calculated using both methods for 
100,000 trees from several different islands. Both methods 
found strongest support for the same tree (Fig. 9; hereafter, 
the ML tree). Traditional tree statistics are presented in 
Table 1 and the figure captions. 

The ML tree is generally similar to the parsimony 
tree (Fig. 1) and also to the estimates provided in Wagner 
(1995a). Trochonema is again placed among Troch- 
onemella spp., whereas other taxa previously assigned to 
the Trochonematidae ( e . g. Eunema, Gyronema , Pro - 
turritella) are derived separately from traditional lophos- 
pirids. This corroborates previous several previous 
hypotheses ( e . g. Ulrich and Scofield, 1897; Wenz, 1938; 
Knight et ai, 1960; Erwin, 1990). 

The ML tree (as the CIS tree before it) derives the 
problematic Lophospira serrulata directly from the con¬ 
temporaneous L. perangulata. That species was previously 
linked by parsimony to the clade of L. rectistriata\ here, 
parsimony links it to the Trochonemella-Trochonema clade. 
L. serrulata shares several homoplasies with other species 
regardless of where it is placed, including open coiling of 
the gerontic whorls, pronounced sutural and umbilical cari- 
nae, and sharp right and left ramp earinae. Notably, these 
features are shared with species known from the same sedi¬ 
ments, such as Trochonema umbilicata, Trochonemella 


Table 1. Properties and statistics of alternative phylogenetic hypotheses. 
(HG, hypothesized gaps in horizon-level sampling; HL, hypothesized tree 
length; L [H I r], likelihood of a hypothesis given the measurable data; 
ML, maximum likelihood; Pars, parsimony; PL, parsimony reconstruction 
of the character matrix; SD, hypothesized stratigraphic debt, with debt 
measured in the same stratigraphic units as “Ranges”). 


Tree 

HL 

L [HL 1 PL] 

SD 

L [SD 1 Ranges] 

L [HG 1 Horizons] 

Pars 

350.50 

<3 x 10' 5 

102 

5.95 x 10' 6 

1.04 x 10' 2U 

ML 

378.37 

0.598 

27 

0.252 

4.96 x 10 43 


notablis, L. helicteres, and Eunema strigillata. These fea¬ 
tures are also homoplastic among these same taxa, which 
makes the assessment of these features as homoplastic on 
L. serrulata all the more plausible. 

Multicharacter homoplasies suggest that characters 
did not evolve independently. However, the homoplastic 
characters appear in different orders among the different 
lineages (/. e. morphotype 11 might be derived from 00 —► 
01 in some cases, 00 —► 10 in others, and directly from 00 
in still other cases). Also, the characters in question all vary 
freely in other parts of the clade. Finally, these characters 
sometimes present the only discernible differences among 
morphospecies. The initial assumptions of character inde¬ 
pendence almost certainly are violated, albeit “fuzzily.” 
Sophisticated likelihood tests examining varying degrees of 
interdependence among characters offer one solution for 
this problem. However, it bears stressing that the likeli¬ 
hood method used here recovers a pattern that violates the 
initial assumptions whereas parsimony does not. 

Regardless of the implications for character coding. 



A Stratigraphic Debt g Stratigraphic Debt 


Fig. 8. A. Estimated sampling intensity (R) given stratigraphic debt distributed among 82 species (with two present in the first interval). Shown are the ana¬ 
lytic solution, the average from 1,000 simulations and the 95%. B. Log-likelihoods of stratigraphic debts (see Fig. 7). 
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Fig. 9. Maximum likelihood estimate of lophospiroid phylogeny. The tree posits 378.4 steps (Cl = 0.293; RI = 0.784) and 25 units of stratigraphic debt, with 
debt units based on the substages marked on the time scale. The tree is considered more likely not simply based on stratigraphic debt (Fig. 10) and continu¬ 
ous stratigraphic data (Table 1), but also on morphologic data (Fig. 8). As in Fig. 1, solid lines show estimated phylogenetic links and dashed lines denote 
ghost lineages or taxa. 


the overall pattern is consistent with a hypothesis of strong 
functional or ecological convergence among these taxa. 
(Lophospiroids lacking these features are known from the 
same strata, so the characters cannot be dismissed simply as 
ecophenotypic variants.) Returning to the particular ques¬ 
tion of Lophospira serrulata , that species retains many ple- 
siomorphic features. Small specimens are very similar to 
L. perangulata , save that the medial lira of the sinus keel is 
serrated and a prominent right ramp carina is present. 
(Note that “sinus keel” is used here instead of selenizone, 
as the latter term is used to denote a morphogenetic artifact 
of a slit. Sinus keels appear long before slits and the fea¬ 
ture represents a separate homology that happens to be in 


the same location as a slit.) If the additional synapomor- 
phies realized on larger specimens are associated with a 
functional complex, then L. serrulata might be separated 
from L. perangulata by fewer evolutionary innovations than 
implied by the character coding. 

Parsimony links Donaldiella bowdeni and other 
high-spired lophospiroids to trochiform species such as 
Lophospira burginensis which differs from the trees pre¬ 
sented by Wagner (1995b). The ML tree differs by (1) 
placing only D. bowdeni and the Silurian D. trilineata in 
that clade, and (2) considering the characters associated 
with the high-spired morphology to be secondarily derived 
rather than primitive. All trees suggest that the earliest 
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lophospiroids were high-spired, corroborating Grabau’s 
(1922) hypothesis about lophospiroid origins. However, 
the ML tree corroborates Ulrich and Scofield’s (1897) 
hypothesis that D. bowdeni was derived from L. oweni 
rather than from older high-spired species. 

LIKELIHOOD RATIO TESTS 

Despite being nearly 10% longer than the most-par¬ 
simonious tree, the ML tree is orders of magnitude more 
likely than the parsimony tree given morphologic data 
alone (Table 1). (Again, the likelihood was calculated 
assuming that character evolution matched the weighting 
scheme of parsimony.) If variable rates of character change 
are used or if patterned homoplasies are included (both of 
which are biologically realistic), then shorter trees become 
even less likely whereas longer ones become more likely. 
The only model of evolution yields most likely lengths 
greater than 370 steps (/. e. the approximate most likely 
length gives a parsimony length of 350.5 steps among inde¬ 
pendently evolving lophospiroid characters): one in which 
there is no homoplasy. However, this hypothesis is falsified 
by the demonstration of incompatibility among lophos¬ 
piroid characters. (In essence, this also falsifies the abduc- 
tive predication of the parsimony proposition “if phyloge¬ 
netic topology, then maximum congruence topology”). 

Likelihood ratio tests can be used to evaluate the 
relative support for competing hypotheses. The test statis¬ 
tic is: 

5 = 2 * (In Li - In L 0 ) 

where L 0 is the likelihood of the null hypothesis (here, the 
parsimony tree) and Li is the likelihood of the test hypothe¬ 
sis (here, the likelihood tree) (Sokal and Rohlf, 1981: 695). 
The results of this test are given in Table 2. Ordinarily 5 is 
evaluated using a chi-square distribution. Goldman (1993) 
cautioned against doing so when evaluating phylogenies 
and suggested using simulations to determine the distribu¬ 
tion of 5 instead. However, simulated 5-distributions differ 
from a chi-square distribution at low 5-values, but not at 
high ones (see Wagner, 1998a). The extremely high 5-val¬ 
ues indicate that both morphologic and stratigraphic data 
provide significantly more support for the ML tree than 
they do for the parsimony tree. 


Table 2. Log-likelihood ratio tests contrasting the parsimony estimate 
with the maximum likelihood estimate. “8” gives the ratio test statistic and 
is evaluated using a Chi-square distribution. 


Stratigraphic Likelihood 

In L[tree 1 PL & Strat. Data] 

Parsimony ML 8 p 

Stratigraphic Debt 
Horizon Data 

-22.45 -1.89 

-496.22 -97.92 

41.11 

796.59 

1.44 x 10' 10 
2.97 x 10" 175 



ML (Continuous) 
ML (Discrete) 

□ Parsimony 
Stratocladi sties 
Cl Sieved 


jsl 


1 2 3 
Number of Incorrect Nodes 


Fig. 10. Performance of maximum likelihood methods. “Continuous” uses 
the method of Huelsenbeck and Rannala (1997) to calculate stratigraphic 
likelihood; “Discrete” uses the method of Foote (1997) to calculate strati¬ 
graphic likelihood. The simulations used the “variable” model of character 
evolution, but the maximum likelihood estimates for morphologic change 
assumed an equiprobable model. Nevertheless, both methods outper¬ 
formed parsimony, stratocladistics, and sieving with confidence intervals. 


THE EFFICACY OF MAXIMUM LIKELIHOOD 
METHODS 

An important consideration when contrasting differ¬ 
ent methods is relative performance under simple cases 
(Felsenstein, 1981, 1984). To this end, I conducted an 
additional set of simulations that contrasted the relative 
efficacy of the two maximum likelihood methods with that 
of parsimony, stratocladistics, and confidence interval siev¬ 
ing. Simulations were conducted using the variable rate 
model of character evolution described above (see Fig. 4). 
Note that the model of character change did not match the 
assumptions of the likelihood estimates of morphologic 
change. In addition, simulated sampling was made more 
variable, with sampling densities and sampling opportuni¬ 
ties varying from interval to interval. 

Both methods perform far better than does parsimo¬ 
ny (Fig. 10), and both also outperform the SCI, CIS, and 
stratocladistics. Huelsenbeck and Rannala’s (1997) method 
does somewhat better than does the debt-likelihood test. 
This is unsurprising because Huelsenbeck and Rannala’s 
(1997) method can distinguish the relative importance of 
particular gaps whereas the former cannot. Interestingly, 
stratocladistics frequently derived the same trees for simu¬ 
lated clades as did maximum likelihood. The stratocladis- 
tic and CIS trees were both very similar to the maximum 
likelihood trees with lophospiroid data. The general con¬ 
gruence among these methods and the superiority of these 
methods over parsimony in simplified simulations both 
suggest that the ML estimate to be a better summary of 
lophospiroid evolution than is the parsimony tree. 
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DISCUSSION 

EVOLUTIONARY IMPLICATIONS OF THE 
PHYLOGENETIC ESTIMATES 

Accurate phylogenies are premises of many tests of 
macroevolutionary hypotheses (Harvey and Pagel, 1991). 
When presented with rival estimates of phylogeny, it is 
important to explore whether the topologic differences 
affect general macroevolutionary interpretations (Donoghue 
and Ackerly, 1996; Wagner, 1997). 

SPECIATION PATTERNS 

Initial phylogenetic estimates suggested that bud¬ 
ding cladogenesis was the most common speciation pattern 
within lophospiroids (Wagner and Erwin, 1995). This 
study also found a positive association between species’ 
durations and numbers of apparent descendant species. 
Few patterns were found consistent with either bifurcating 
cladogenesis or anagenesis, both of explicitly predict that 
ancestral and derived species did not co-occur. (Note that a 
“species” here is a unique combination of character states.) 
Budding cladogenesis does not demonstrate any particular 
speciation process. However, it is consistent with processes 
such as peripheral isolation (Mayr, 1963), shifting-balances 
(Wright, 1931; but see Provine, 1986), and punctuated 
equilibrium (Gould, 1982) and inconsistent with processes 
such as selective divergence (Darwin, 1859) and vicariance 
(Brooks and McLennan, 1991). 

Both the parsimony and ML trees corroborate the 
original conclusions about speciation patterns. There are 
18 taxa with plesiomorphic sister species on the parsimony 
tree. In 15 of those cases, the plesiomorphic taxon co¬ 
occurs with the derived sister taxon. There are 49 species 
with plesiomorphic sister species on the ML tree, 36 of 
which co-occur with their putative ancestor. On both trees, 
significantly more cases are consistent with budding clado¬ 
genesis than with anagenesis or bifurcating cladogenesis 
(parsimony: p- 1.32 x 10’ 3 ; ML: p = 1.42 x 10" 4 , based on 
a standard binomial test assuming a 50:50 distribution). 
However, one major difference exists. The parsimony tree 
posits no species with multiple descendants whereas the 
ML tree posits several species with two or more descen¬ 
dants. Moreover, there is a strong association between 
duration and the apparent number of descendants on the 
ML tree (Kendall’s x= 0.391; p = 2.44 x 10' 7 ; Sokal and 
Rohlf, 1981: 429). Models predicting budding cladogenic 
patterns are generally consistent with the pattern shown by 
the ML tree. However, it is not clear that any speciation 
models predict the pattern shown by the parsimony tree, 
where species such as Lophospira perangulata produce a 
single daughter taxon early and then persist many millions 
of years with no additional daughter taxa. 


TRENDS 

Wagner (1996) documented active trends for several 
shell features among Ordovician-Silurian gastropods, 
including lophospiroids. One of these trends, the reduction 
of the sinus, is clearly observable among lophospiroids. 
Early species such as Lophospira perangulata have wide, 
deep sinuses that curve continuously back toward the sinus 
keel. Similar sinuses are retained among species classified 
as Donaldiella and in some Lophospira subcladcs ( e . g. the 
L. burginensis clade). However, reduction of the sinus hap¬ 
pens in parallel among derived TrochonemeUa (/. e. the 
clade including T. montrealensis and derived Trochonema 
(/. e. the clade including T. umbilicata). Note that the com¬ 
mon ancestor of those two clades (/. e. species similar to 
TrochonemeUa knoxvillensis and T. trochonemoides ), 
retained a L. perangulata -like sinus. A similar pattern of 
sinus reduction is observed among species in the clade 
including L. tnilleri). It should be noted that the trend as 
described at this level can be inferred from the parsimony 
and ML trees. 

Additional parallelisms concerning the sinus exist 
on both estimated trees. In the case of derived Trochonema 
and the clade including Lophospira centralis , sinuses lose 
most of their curvature and retreat nearly straight to the 
sinus keel. A parallel reduction of the sinus keel is 
observed within both clades. Lophospiroids primitively 
possess a trilineate sinus keel, with a strong, sharp, medial 
lira bordered by two sharp and somewhat weaker peripheral 
lira. Broader analyses [Wagner, in press (a)] imply that the 
peripheral lirae of lophospiroids are homologous with the 
peripheral lira of other early Murchisoniinae (e. g. 
Hormotoma and Eotomaria) whereas the medial lira is a 
synapomorphy of early lophospiroids. The peripheral lirae 
are lost on Trochonema and twice lost within the L. cen¬ 
tralis clade (/. e. Eunema and Gyronema spp.). The periph¬ 
eral lirae also are lost in Proturritella and the 
Kiviasukkaan-Loxoplocus clade. Thus, the single peripher¬ 
al keel of species previously assigned to the 
Trochonematoidea (e. g. Knight et al ., 1960) appears to be 
the medial lira of primitive lophospiroids. The reduction of 
the sinus keel precedes the reduction of the sinus among 
Trochonema spp., because species such as T. bellula and T. 
pandori possess only the medial lira but retain curved and 
fairly wide sinuses. Conversely, the reduction of the sinus 
precedes the reduction of the sinus keel among the L. cen¬ 
tralis clade, because species such as L. centralis and L. 
helicteres possess shallow, straight sinuses but retain a tri¬ 
lineate sinus keel. Similar mosaic patterns of sinus and 
sinus-keel reduction are observed within other gastropod 
clades [Wagner, in press (a)]. Thus, the characters clearly 
are independent homologues in that a change in one does 
not necessitate an immediate change in the others. 
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However, the overall pattern suggests some sort of biologi¬ 
cal (e. g. functional) association. 

Bandel and Geldmacher (1996) have stated that 
species such as Trochonema and Eunema could not be 
derived from ancestors with sinuses and/or slits and thus 
could not be lophospiroids. However, early Trochonema 
spp. possess a lophospiroid sinus that is reduced or entirely 
lost in derived species. Eunema spp. retain a derived sinus 
observed in species previously assigned to Lophospira 
(e. g. “L.” centralis). Sinus reduction obviously was com¬ 
mon among Ordovician and Silurian lophospiroids. A slit 
is an uncommon feature among lophospiroids, but likely 
evolved at least three times within the clade. Derived 
Trochonemella spp. such as T. notablis are one such group 
of taxa, but primitive species such as T. ktioxvillensis and T. 
trochonemoides clearly did not possess slits. Thus, whether 


Trochonema or Eunema could have slit-bearing ancestors is 
irrelevant, as neither ML nor parsimony trees posit slit¬ 
bearing ancestors for either taxon. 

Parsimony and ML trees suggest very different pat¬ 
terns of ornament evolution. The increasing number of 
ornate gastropod taxa over the Phanerozoic is a well-docu¬ 
mented trend (e. g. Vermeij, 1977). Ornament is uncom¬ 
mon among lophospiroids until the Late Ordovician. 
Conversely, most Silurian lophospiroids had at least some 
ornament. The parsimony tree suggests that ornament 
arose four times and was subsequently lost four times (Fig. 
11). (The reconstructions on the parsimony tree are unam¬ 
biguous and thus not affected by optimization assump¬ 
tions.) Parsimony also suggests that ornate taxa were sub¬ 
stantially more diverse in the Middle Ordovician than sam¬ 
pling alone would imply. Hypothesized ornate range exten- 



Fig. 11. Patterns of ornament evolution as posited by parsimony. Black branches denote primitively inornate species, gray branches denote secondarily inor¬ 
nate species, and white branches denote ornate species. Note that the tree posits a much earlier divergence of ornate species than observed in the fossil 
record, resulting in hypotheses of many ornate ghost taxa (/. e. white dashed lines). 
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sions (Fig. 11; dashed white lines) represent approximately 
one-quarter of the inferred diversity during the Llanvirn. 
However, only one ornate species is known before the 
Middle Caradoc. Parsimony also implies that several inor¬ 
nate Ordovician species are secondarily inornate. 
Suspiciously, ornate relatives always appear after the “sec¬ 
ondarily” inornate taxa. 

The ML tree suggests that ornament arose nine 
times (Fig. 12), which implies a frequency of change of 
0.098 per branch. This is over twice the typical frequency 
(0.040 per branch). However, ornament apparently is local¬ 
ly conservative, as there are no apparent reversals in sub¬ 
codes of up to seven ornate species. This pattern suggests 
a driven trend (/. e. the biased production of a particular 
morphologic type; McShea, 1994). Despite the low sample 


size, a G-test on exact binomial probabilities rejects any 
hypothesis of gains:losses ratios lower than 77:23. Vermeij 
(1977) hypothesized that increasing predation pressures 
drove a trend towards greater proportions of ornate gas¬ 
tropods. Taxic data suggest that this trend was well under¬ 
way by the Devonian (Signor and Brett, 1984). 
Lophospiroids suggest that the trend began by the Late 
Ordovician, implying that either (1) predation pressures 
began earlier than previously suspected, or (2) the trend 
was driven (at least initially) by factors other than preda¬ 
tion. 

There are reasons why parsimony would have prob¬ 
lems reconstructing Fig. 12. First, simulation studies indi¬ 
cate that frequencies of change half that of ornament char¬ 
acters will induce errors in parsimony (Mooers et al ., 



Fig. 12. Patterns of ornament evolution as posited by likelihood. Black branches denote primitively inornate species and white branches denote ornate 
species. No secondarily inornate species are posited. The evolution of ornate species per substage here matches patterns observed in the fossil record and 
suggests a driven trend (sensu McShea, 1994). 































































16 


AMER. MALAC. BULL. 15(1) (1999) 


1995). Simulations further suggest that parsimony is even 
less successful when transition patterns are biased (Kuhner 
and Felsenstein, 1994; Lamboy, 1994). Second, it is 
impossible to code ornament characters in a way that is 
neutral with regard to all possible patterns of evolution. 
Some species possess ornament throughout the shell, but 
species such as Ruedemannia lirata and Arjamannia 
thraivensis possess ornament only on the left ramp. In 
some non-lophospiroids, ornament occurs only on the right 
ramp. Thus, the presence and absence of ornament must be 
scored separately for both ramps. The result is that com¬ 
pletely ornate species share two synapomorphies instead of 
one; if evolution goes from completely inornate to com¬ 
pletely ornate (e. g. Lophospira milleri to Proturritella 
bicarinata ), then parsimony counts two steps instead of 
one. A third and related reason is that similar characters 
were coded as synapomorphic (sensu Sober, 1988), if only 
at basic levels. For example, a species with coarse, broad¬ 
ly-spaced ornament on the left ramp shared the presence of 
left ramp ornament (character 87) with a species with fine, 
densely-packed ornament. Additional characters ( e . g . 
characters 88 and 89) then were used to distinguish the two. 
This assumes that it is easier to modify a character than it is 
to evolve it, which many workers regard as the lesser 
assumption ( e . g. Hennig, 1966). 

The trend toward increasing numbers of ornate 
species highlights a potential strength of the maximum like¬ 
lihood method over parsimony. Parsimony is ill-equipped 
to deal with biased transitions and character transitions that 
cannot be neatly coded. As such, the pattern of evolution 
shown in Fig. 12 predicts that parsimony results such as 
shown in Fig. 11 should be probable. However, the maxi¬ 
mum likelihood analysis (with the help of stratigraphic 
data) was able to recover a pattern of evolution that was at 
odds with the initial assumptions. 

The effect of initial assumptions about character 
change has additional implications. The ML tree suggests a 
trend towards increasing numbers of ornate species induced 
by a tendency for inornate ancestors to yield ornate descen¬ 
dants. The parsimony tree suggests a trend driven by a very 
speciose ornate clade. Several workers ( e . g. Wright, 1967; 
Eldredge and Gould, 1972; Stanley, 1975; Vrba and 
Eldredge, 1984) hold that differential speciation might 
induce trends and that morphologic evolution might be ran¬ 
dom with respect to those trends (/. e . “Wright’s Rule”; 
Gould, 1977). By minimizing the number of times charac¬ 
ters arose and by maximizing the number of taxa linked by 
those characters, parsimony might be biased towards results 
consistent with “species selection” rather than with driven 
trends. Lophospiroids offer a possible example of this bias. 

RECOMMENDED TAXONOMIC REVISIONS 

Phylogenetic frameworks serve as a basis for taxo¬ 


nomic classification, which in turn can convey historical 
patterns of biodiversity to non-specialists. Conveying 
diversity patterns taxonomically is not simple. Some work¬ 
ers ( e . g. Eldredge and Novacek, 1985; Smith and 
Patterson, 1988) argue that only monophyletic supraspecif- 
ic taxa can reflect historical diversity patterns. If so, then 
phylogenetic systematics (Hennig, 1966) is appropriate. 
However, other workers use higher taxa as proxies for 
species in diversity studies because higher taxa are less 
prone to sampling biases than are species ( e . g. Raup, 1979; 
Sepkoski, 1979, 1993; Raup and Boyajian, 1988). 
Paraphyletic and monophyletic taxa both are excellent 
proxies for extant species diversities (Williams and Gaston, 
1994; Andersen, 1995; Balmford et al ., 1996). Both simu¬ 
lation (Sepkoski and Kendrick, 1993; Maley et al ., 1997) 
and empirical studies (Wagner, 1995c) indicate that taxo¬ 
nomic schemes including paraphyletic taxa portray species- 
level diversity patterns as well as or better than do those 
using monophyletic taxa alone. Including paraphyla 
becomes especially important when rates of cladogenesis 
and/or extinction or sampling intensities vary over time 
(Wagner, 1995c). Thus, a taxonomic scheme designed to 
reflect historical diversity patterns should include para¬ 
phyletic taxa lest it be biased against depicting particular 
patterns. 

Evolutionary systematics (Simpson, 1961) allows 
for paraphyletic taxa, but evolutionary systematics has no 
explicit, repeatable criteria for delimiting such taxa (Hull, 
1967; note that the same actually is true of Hennigian sys¬ 
tematics; Tassy, 1988). Phylogenies are nothing more than 
artifacts of speciation and extinction rates (Sepkoski, 1987). 
Such an evolutionary parameter can be used to delimit 
“real” groups in a statistical sense. For example, quantita¬ 
tive analyses of net rates of cladogenesis in primates 
(Purvis et al ., 1995) and rates of morphologic change in 
rostroconch molluscs (Wagner, 1997) separate paraphyletic 
groups from monophyletic ones. Unfortunately, demon¬ 
strating that sets of species belong to statistically distin¬ 
guishable groups requires large sample sizes, much larger 
than is available for defining genera or probably even fami¬ 
lies. Estabrook (1986) proposed integrating phenetic dis¬ 
similarity with phylogenetic topology to delimit higher 
taxa. This has the advantage of being fully repeatable once 
criteria are made explicit. However, if rates of morphologic 
evolution vary over time or across a clade, then more higher 
taxa will unite very different numbers of species and thus 
will not portray historical biodiversity accurately. 

Although based on the ML tree, the generic defini¬ 
tions presented here would be nearly identical if derived 
from the parsimony tree. The lack of explicit criteria for 
delimiting genera does not interfere with a generic revision 
that better reflects historical diversity patterns than current 
generic schemes do (Fig. 13). For example, Lophospira 
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Fig. 13. Proposed taxonomic revisions. denotes the type species of each genus. “Proturritella ” is retained because additional species not available for 
this analysis might belong there. Formal designation of “?New Genus” is withheld pending additional analyses, because it is possible that those species are 
instead affiliated with Donaldiella filosa, the type species of Donaldiella. If so, then the species here placed in Donaldiella should be returned to 
Pagodospira. See Appendix 4 for necessary rediagnoses and redefinitions. 


bellicarinata could be assigned to ether Arjamannia or 
Loxoplocus ; in either case, one genus appears in the Late 
Caradoc whereas the other appears in the Ashgill. 
However, retaining all Arjamannia and Loxoplocus spp. in 
Lophospira results in a taxonomy that ignores the deriva¬ 
tion of a diverse clade whose sister taxa (/. e. other deriva¬ 
tives of L. milleri such as Ruedetnannia and the Eunema- 
Gyronema clade) are separated into one or more genera. 

In some cases, the revisions simply update lophos- 
piroid taxonomy. Many of the species analyzed here were 
originally classified before Lophospira was partitioned into 
additional genera. Initial diagnoses of new genera were not 
accompanied by exhaustive revisions of the entire clade, 


resulting in a highly polyphyletic Lophospira. This is 
remedied to make Lophospira a paraphyletic group. Other 
cases involve redefining and rediagnosing obviously poly¬ 
phyletic genera. For example, Eunema previously has been 
diagnosed by high translation (sensu Raup, 1966), the 
absence of a sinus and a medial keel. A sinus actually is 
present in the type species, albeit a shallow V-shaped one 
that is shared with some species formerly placed in 
Lophospira (all now reclassified as Eunema). In the redefi¬ 
nition presented here, only the type species lacks a trilin- 
eate sinus keel. However, the loss of a peripheral lira repre¬ 
sents a derived condition that happened multiple times and 
therefore cannot unequivocally diagnose genera. 
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Paupospira 


“ Lophospira ” "Lophospira" “ Lophospira ” “Donaldiella" “Lophospira” “Lophospira" “Schizolopha” 

burginensis tennuistriata oweni bowdeni sumnerensis tropidophora moorei 



Fig. 14. Cladogram of Paupospira highlighting the important morphologic features of the clade. White branches denote zero-length branches (/. e. species 
identical to hypothesized ancestral morphotypes). 


Patterson and Smith (1987) suggested that some 
diversity patterns inferred from taxic data might be artifacts 
of paraphyly. However, if a paraphyletic taxon’s disappear¬ 
ance is not due to its sole species anagenetically transform¬ 
ing into a species placed in a different higher taxon, then its 
disappearance necessarily reflects the extinction of at least 
one species. Because paraphyletic taxa necessarily are 
older than their monophyletic “descendants,” paraphyla 
often have higher standing diversities than their mono¬ 
phyletic “daughters” (Uhen, 1996) and thus require greater 
numbers of species extinctions to be terminated. For exam¬ 
ple, Trochonemella is paraphyletic relative to Trochonema , 
which in turn is paraphyletic relative to Globonema. 
However, the extinctions of Trochonemella and 
Trochonema in the Ashgill represents the loss of eight 
species (as well as the disappearance of particular morpho¬ 
types). Globonema represents two survivors, so the 2:1 
casualties.survivors ratio implied by genera underestimates 
the 4:1 ratio implied by species. Alternative classifications 
fare worse. Any monophyletic definition including all 
members of the clade fails to register the extinction. 


Conversely, reduced monophyletic definitions recognizing 
only the clades terminated in the Ashgill fail to recognize 
that the group diverged during the Middle Ordovician. A 
new monophyletic genus, Paupospira , also illustrates how 
strictly monophyletic higher taxa can obfuscate historical 
patterns. That genus is nearly eliminated by the end- 
Ordovician extinction, but a sole survivor is known from 
the Llandovery. The genus’ demise somewhat misleadingly 
now represents a “background” extinction (Jablonski, 
1986). 

CONCLUSIONS 

Although malacologists have questioned the utility 
of shell characters in phylogenetic analyses, shell character 
data for Ordovician-Silurian lophospiroids show compati¬ 
bility similar to that expected if those data do retain hierar¬ 
chical (phylogenetic) signal. Nevertheless, truly hierarchi¬ 
cal signal among the clade as a whole and geologically 
older taxa in particular decreases with the addition of geo¬ 
logically younger taxa. Phylogenetic methods incorporat- 
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ing stratigraphic data outperform parsimony for simulated 
clades with similar fossil records and the same number of 
characters and character states as lophospiroids. 
Correspondingly, the parsimony estimate of lophospiroid 
phylogeny posits many more gaps than plausible given the 
quality of the lophospiroid record. 

Likelihood tests support a tree that is 30 steps 
longer but nevertheless is more likely based on both mor¬ 
phologic and stratigraphic data. Both trees suggest that 
budding cladogenesis was the dominant speciation pattern, 
and both portray trends toward increasing numbers of 
ornate species. However, the likelihood tree depicts a 
strong association between duration and proclivity at the 
species level whereas the parsimony tree does not. Also, 
the likelihood tree suggests that the ornament trend was dri¬ 
ven by the acquisition of ornament being far more common 
that its loss. Parsimony suggests that ornate clades were 
simply more diverse that were inornate clades. At least in 
this example, parsimony appears biased toward results sug¬ 
gesting patterns consistent with “species-selection.” 
Lophospiroids highlight the utility of gastropods as a model 
group for testing macroevolutionary hypotheses. However, 
lophospiroids also highlight the importance of developing 
phylogenetic methods in whose results we can be more 
confident. 
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APPENDIX 1 

CHARACTERS AND CHARACTER STATES 


Character state codings follow Wagner [in press (a)]. Lophospirioids do 
not possess all of the states available in some multistate characters devel¬ 
oped for that analysis; therefore, in the following list, the plesiomorphic 
state for lophophorids is not coded “1” in all cases, and some character 
states appear to be missing. The character coding has not been adjusted 
herein, so that lophospiroids can be easily compared with other gas¬ 
tropods. 

1. Sinus (broad cleft in shell, culminating at presumed location of anus): 

(1) absent; (2) present. Weight 1. 

2. Sinus depth [as described by general angle of retreat by sinus; continu¬ 

ous character, ranging from almost 0° to 40°; really a shape metric, 
as true depth is a product of angle of retreat, sinus width (character 
3), and curvature of retreat (character 4)]: (1) >0°; (2) -10°; (3) 
-20°; (4) -30°; (5) -40°. Ordered; weight 0.25. 

3. Sinus width [as described by angle passing from center of aperture 

through top and bottom of sinus; continuous character, with states 
broken into just above/below sinus apex (1), halfway between 
top/bottom of ramp and sinus apex (2), coinciding with top/bottom 
of ramp (3), and beyond top/bottom of ramp (4); coding the position 
of the onset of the sinus relative to the top/bottom of the ramps rather 
than the absolute width accommodates the fact that changes in ramp 
lengths and orientations also will change sinus width as an artifact; 
this coding scheme focuses only on differences clearly attributable to 
the sinus itself]: (1) just above SK; (2) between SK and RR; (3) at 
RR. Ordered; weight 0.5. 

4. Sinus shape [describing general trend of sinus; nearly straight up to 

sinus apex (/. e. V-shaped) or curving back continuously; this 
describes the shape of the curve whereas sinus depth (character 2) 
describes the general angle of retreat]: (2) straight; (3) continuous. 
Weight 0.5. 

5. Crenulated aperture [in which apertural margin generates “zig-zag” 

growth lines]: (1) absent; (2) present. Weight 1. 

6. Sigma-shaped aperture: (1) absent; (2) intersection between base and 

alveozone curved (yielding reverse sigma) instead of planar or angu¬ 
lar. Weight 1. 

7. Prominence of growth lines (GL) [evaluated only by comparison with 

both different specimens known from same preservational conditions 
and similar specimens from different preservational conditions]: (2) 
weak; (3) fine sharp; (4) strong. Ordered; weight 0.5. 

8. Imbricated GLs [alternating growth lines appreciably stronger than oth¬ 

ers; considered present only if imbrications obviously are patterned, 
as random imbrications can occur due to shell repair or on gerontic 
whorls of some species; growth line prominence is coded based on 
weaker growth lines]: (1) absent; (2) weak; (3) moderate; (4) strong. 
Unordered; weight 1. 

9. Sinus keel (SK) width [originally thought to be morphologic artifact of 

a slit or a “notch,” however, appearing on slitless specimens and the 
seemingly mythical “notch” has never been documented; continuous 
character, where 360( equals circumference of aperture]: (1) ~5°; (2) 
-10°; (3) -15°; (4) -20°; (5) -25°; (6) -30°. Ordered; weight 0.2. 

10. Peripheral lirae (PL) [two lirae on bilineate peripheral bands (or outer 

two on trilineate bands of lophospiroids); in some cases, might be 
produced by placing a slit within a medial lira (see below); however, 
not occuring on all species with slits and many species with periph¬ 
eral lira clearly do not have slits]: (1) absent; (2) present. Weight 1. 

11. PL type: (1) threads with round profile ( e . g. Ruedemannia)’, (2) 
threads with sharp profile ( e . g. Ectomaria or Lophospira). Weight 1. 

12. PL strength: (2) weak (clearly visible but casting little relief); (3) 
moderate (casting noticeable relief); (4) strong (clearly projecting 
from shell). Ordered; weight 0.5. 


13. Medial lira (ML) [= “notch keel” or “sinus keel”; present with periph¬ 

eral lirae (character 10) only if there is ontogenetic change (character 
25) or if peripheral band is trilineate; carina at the apex of the sinus]: 

(1) absent; (2) present. Weight 1. 

14. ML type [as seen in profile]: (1) round, of equal “height” and width 

( e . g. Ruedemannia)\ (2) sharp, of equal “height” and width (e. g. 
Lophospira). Weight 1. 

15. ML strength: (1) extremely weak (barely visible, casting almost no 
relief); (2) weak (clearly visible but casting little relief); (3) moder¬ 
ate (casting noticeable relief); (4) strong (clearly projecting from 
shell). Ordered; weight 0.33. 

16. Imbricated ML; (1) consistent; (2) flaring periodically ( e . g. 
Lophospira serrulata). 

17. SK prominence [not equal to strength, as prominence refers to whole 

structure projecting from rest of shell; sometimes coincides with 
channel underneath peripheral band; however, species with channels 
and non-prominent peripheral bands and species with prominent 
peripheral bands and no obvious channel are both known]: (1) 
peripheral band not altering profile of whorl; (2) whole peripheral 
band projecting slightly, creating slight ridges between shell and 
peripheral band; (3) whole peripheral band projecting strongly, creat¬ 
ing strong ridges between shell and peripheral band. Ordered; 
weight 0.5. 

18. Slit: (1) absent; (2) present. Weight 1. 

19. Slit depth: (1) shallow (extending <10° behind aperture, measured 
looking down coiling axis); (2) deep (extending -20° behind aper¬ 
ture). Weight 1. 

20. Maintenance of slit: (1) periodically erased or reduced, with shell 
material deposited within slit less frequently than on rest of shell, 
resulting in slit depth (and sometimes presence/absence) varying 
over time ( e. g. Trochonemella\ recognized in part by greater dis¬ 
tance between lunulae than observed between growth lines); (2) con¬ 
tinuously maintained {e. g. ,, Schizolopha’ > moorei). Weight 1. 

21. Lunulae shape [lunulae here are growth lines within sinus keel]: (1) 

concentric (shallow U-shape); (4) V-shaped ( e . g. Lophospira cen- 
tralis and relatives). Weight 1. 

22. Lunulae strength: (1) weaker than GLs; (2) same as GLs; (3) 
stronger than GLs (/. e. imbricated; e. g. Trochonemella spp.). 
Unordered; weight 1. 

23. Imbricated lunulae type: (1) deep and obtuse; (2) serrated. Weight 1. 

24. Ontogenetic change in imbricated lunulae: (1) juvenile whorls only; 

(2) throughout. Weight 1. 

25. Reduction to monolineate SK over ontogeny: (1) absent; (2) monolin- 

eate on final whorls. Weight 1. 

26. Peripheral band attitude: (1) projecting straight from shell, as seen on 

most species; (2) curving adapically, e. g. Arjamannia and relatives. 
Weight 1. 

27. Midwhorl (MW) channel [groove underneath medial lira; often but 
not always associated with prominent band, and so coded separate¬ 
ly]: (1) absent; (2) present. Weight 1. 

28. MW channel strength: (1) weak; (2) strong. Weight 1. 

29. SK position [described relative to aperture centroid, with 90( indicat¬ 

ing that a plane is passing through centroid and that SK is perpendic¬ 
ular to coiling axis; 0( indicating that plane is parallel to coiling 
axis]: (2) -100°; (3) -90°; (4) -80°; (5) -70°; (6) -60°; (7) -50°. 
Ordered; weight 0.2. 

30. Ramp shape symmetry [on primitive bilaterally symmetrical gas¬ 
tropods, right and left ramp shapes are symmetrical; derived species 
are asymmetrical]: (1) right rounder; (2) symmetrical; (3) left 
rounder. Unordered; weight 1. 
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31. Right ramp (RR) shape: (1) globular; (2) convex; (3) flat; (4) slightly 

concave; (5) concave. Ordered; weight 0.125 (reflecting both con¬ 
tinuous nature and de-weighting for asymmetry). 

32. RR and LR lengths [on primitive bilaterally symmetrical gastropods, 

right and left ramp lengths are equal ; lophospirds with longer right 
and longer left ramps are known]: (1) longer RR; (2) equal lengths; 

(3) longer LR. Unordered; weight 1. 

33. RR length [describing angle from “top” of ramp to sinus apex, based 

on triangle passing through aperture centroid]: (3) -50°; (4) -60°; 
(5) ~70°; (6) -80°; (7) -90°. Ordered; weight 0.25. 

34. LR length [see character 33]: (3) -50°; (4) -60°; (5) -70°; (6) -80°; 

(7)-90°. Ordered; weight 0.25. 

35. RR:LR projection [ramp projection describing angle of “rise” of ramp 

from sinus keel, based on plane passing through SK and aperture 
centroid]: (1) RR projection higher; (2) LR and RR projections 
equal. Weight 1. 

36. RR projection [see character 35]: (5) -60°; (6) -70°; (7) ~80°(. 
Ordered; weight 0.25. 

37. LR projection [see character 35]: (3) -40°; (4) -50°; (5) -60°; (6) 
-70°; (7) -80°. Ordered; weight 0.125. 

38. Sutural canna (SC): (1) absent; (2) present. Weight 1. 

39. SC strength: (1) weak (creating only a weak profile); (2) moderate 
(partly filling up the suture). 

40. RR carina (RRC) [strong carina usually located at top of sinus]: (1) 
absent; (2) present. Weight 1. 

41. RRC strength: (1) weak (creating no profile); (2) moderate (roughly 
equal to weak-to-moderate peripheral band); (3) strong (roughly 
equal to strong peripheral band). Ordered; weight 0.5. 

42. RRC type: (2) thin local thickening; (3) round profile; (4) sharp pro¬ 

file. Unordered; weight 1. 

43. Ontogenetic change in RRC strength: (1) absent; (2) becoming weak¬ 

er on adult whorls. Weight 1. 

44. Channel beneath RC: (1) absent; (2) present. Weight 1. 

45. RRC location: (1) -75° toward suture from SK; (2) -45° toward 
suture from SK. Weight 1. 

46. RRC attitude: (1) carina projecting perpendicularly to RR; (2) carina 

curving abapically. Weight 1. 

47. RRC: (1) plain thread; (2; serrated. Weight 1. 

48. Shape of shell at top of right ramp [usually at suture]: (1) oblique; (2) 

acute, with channel underneath suture. Weight 1. 

49. LR shape [see character 31 ]: (1) globular; (2) very convex; (3) slight¬ 

ly convex; (4) flat. Ordered; weight 0.167. 

50. Swelling at intersection of left ramp and base [might be primitively 
homologous with left ramp carina (character 51), and present on ear¬ 
liest gastropods]: (1) absent; (2) present. Weight 1. 

51. LR carina (LRC): (1) absent; (2) present. Weight 1. 

52. LRC type: (1) thick contusion, but creating more distinct profile than 

simple swelling; (2) sharp profile. Weight 1. 

53. LRC strength [see character 41]: (1) weak; (2) moderate (~ sinus); (3) 

strong. Ordered; weight 1. 

54. LRC [see character 47; species with serrated right carinae but plain 
left carina exist, indicating that the two evolve independently]: (1) 
plain thread; (2) serrated. Weight 1. 

55. Second LRC: (1) one carina; (2) two carinae ( e . g. Lophospira 
quadrisulcata). Weight 1. 

56. Channel beneath LRC [see character 44]: (1) absent; (2) present. 
Weight 1. 

57. Columella thickness: (1) no thicker than rest of shell; (2) slightly 
thicker than rest of shell; (3) much thicker (partly filling umbilicus); 

(4) extremely thick (filling umbilicus). Ordered; weight 0.33. 

58. Umbilical carina (UC) [carina at base of shell, circling umbilicus]: 
(1) absent; (2) present. Weight 1. 

59. UC type: (1) thick, dull protrusion; (2) sharp extension accomodating 


a channel; (3) lirum with sharp profile. Unordered; weight 1. 

60. UC strength [see characters 41 and 53]: (1) very weak; (2) weak; (3) 

strong. Ordered; weight 0.5. 

61. Ontogenetic change in UC strength: (1) constant; (2) becoming weak¬ 

er on adult whorls. Weight 1. 

62. UC location [with larger angle indicating UC closer to coiling axis]: 
(1) -120° below SK; (2) -90° below SK. Weight 1. 

63. Angle at base of columella [narrower angle indicates sharper, more 
siphonate base of shell]: (3) -60°; (4) -75°; (5) -90°; (6) -105°; (7) 
-120°. Ordered; weight 0.25. 

64. Shape of columella on inner margin [can differ from shape on outer 
margin, so the two coded separately]: (1) arching like half-circle; (2) 
arching in obtuse curve; (3) trending toward straight; (4) curving 
slightly into the aperture. Ordered; weight 0.5. 

65. Outer margin shape: (1) more obtuse than inner margin; (2) same as 

inner margin; (3) more acute than inner margin. Weight 1. 

66. Ontogenetic change in margin shape: (1) none; (2) becoming rounder 

over ontogeny. Weight 1. 

67. Columella attitude [describing main trend of columella relative to 
coiling axis]: (1)0° (/. e. perpendicular to coiling axis); (2) 15°; (3) 
30 c ; (4) 45°; (5) 60°. Ordered; weight 0.25. 

68. Columella lira [carina in middle of columella, visible in umbilicus]: 
(1) absent; (2) present. Weight 1. 

69. Parietal inductura thickness [silicification can occur differently among 

different shell layers, which can affect characters such as this; there¬ 
fore, “relative” states were coded based on comparisons among taxa 
found from the same beds, with comparisons among conspecifics 
from different beds used to establish the final character code]: (1) 
absent; (2) thinner than rest of shell; (3) thickness same as shell. 
Unordered; weight 1. 

70. Columella reflected around coiling axis: (1) absent; (2) present. 
Weight 1. 

71. Whole aperture inclined (tangential): (1) absent; (2) present. Weight 

1. 

72. Degree of whole aperture inclination (InAn): (1) -10°; (2) -20°; (3) 

-30°; (4) -40°; (5) -50°; (6) -60°. Ordered; weight 0.2. 

73. Left side of aperture only inclined [or inclined at different angle than 

right side]: (1) absent; (2) present. Weight 1. 

74. Left InAn: (1)-10°; (2)-20°; (3)-30°. Ordered; weight 0.5. 

75. Right side of aperture only inclined [or inclined at different angle than 

left side]: (1) absent; (2) present. Weight 1. 

76. Right InAn: (1) -10°; (2) -20°; (3) -30°. Ordered; weight 0.5. 

77. Anterior projection of aperture [best observed from base, because 
growth lines will project forward instead of radially; it is important 
that they do this at the onset, as tangential aperture will cause growth 
lines on the periphery to slope forward as well]: (1) absent; (2) pre¬ 
sent. Weight 1. 

78. Degree of anterior projection [relative to plane passing radially 
through coiling axis]: (1) 10°; (2) 20°; (3) 30°; (4) 40°; (5) 50°. 
Ordered; weight 0.25. 

79. Aperture expansion [expansion of shell in radians as a “tube”]: (2) 
0.05-0.10; (3) 0.10-0.15; (4) 0.15-0.20. Ordered; weight 0.5. 

80. Curvature about coiling axis [in radians]: (3) 0.65 < K < 0.75; (4) 
0.75 < K < 0.85; (5) 0.85 < K < 0.95. Ordered; weight 0.5. 

81. Translation (T) [vector of downwards growth (in radians)]: (5) 0.16 < 

T < 0.47 (low dextral); (6) 0.47 < T < 0.79 (moderate dextral); (7) 
0.79 < T < 1.10 (high dextral); (8) 1.10 < T < 1.41 (very high dex¬ 
tral). Ordered; weight 0.33. 

82. Continuous ontogenetic change in T: (2) isometric; (3) continuously 

increasing, resulting in ever-increasing apical angle. Weight 1. 

83. Early ontogenetic change in T: (1) early decrease of T ( e . g. 
Trochonema, Trochonemella)\ (2) isometric; (3) early increase in T. 
Unordered; weight 1. 
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84. Late ontogenetic change in T [often resulting in open-coiled “cork- 
screw”-like shells, e. g. Lophospira helicteres, Eunema strigillata]: 
(2) isometric; (3) late increase in T. Weight 1. 

85. Punctuated late ontogenetic change in T: (1) occuring over 1-2 
whorls ( e . g. Trochonema, Trochonemella)', (2) punctuated, occuring 
over less than one revolution. Weight 1. 

86. Magnitude of late ontogenetic change in T: (1) slight, changing 
suture point; (2) major, resulting in open coiling. Weight 1. 

87. Ornament on left side of aperture: (1) absent; (2) present on left ramp 

and base; (3) present on left ramp only. Unordered; weight 0.5. 

88. LR ornament density [based on average angular distance between 
threads]: (1) 1 per 20°; (2) 1 per 10°; (3) 1 per 5°; (4) 1 per 1°. 
Ordered; weight 0.5. 

89. LR ornament type: (1) thin local thickenings, little stronger than 
growth lines and with no profile ( e . g. Ruedemannia)', (2) lirae with 
weak profile ( e . g. Arjamannia)', (3) thick lirae with strong profile ( e . 


g. Longstajfia). Weight 1. 

90. RR ornament [because species with left ramp ornament sometimes 
lack right ramp ornament, these characters coded separately; com¬ 
pletely inornate to completely ornate represents one step] (1) 
absent; (2) present. Weight 0.5. 

91. RR ornament density [see character 88]: (1)1 per 20°; (2) 1 per 10°; 

(3) 1 per 5°; (4) 1 per 1°. Ordered; weight 0.5. 

92. RR ornament strength: (1) thin threads; (2) weak lirae; (3) strong 
lirae. Ordered; weight 0.5. 

93. Ornament as changes in aperture shape: (1) simple local thickenings 

of shell; (2) representing local flarings of aperture. Weight 1 

94. Ontogenetic change in RR ornament: (1) constant; (2) weaker on 
adult whorls. Weight 1. 

95. Size: (1) <10 nun 3 (micro-mollusk); (2) small (>10 to <10 2 mm 3 ); (3) 

moderate (>10 2 to <10 3 mm 3 ); (4) large (> 10 3 to <10 4 mm 3 ). 
Ordered; weight 0.33. 


APPENDIX 2 
CHARACTER MATRIX 

Character data for lophospiroid species. See Appendix 1 for characters and character descriptions. [?, characters that could not be observed; characters 
that do not pertain to a species ( e . g. ornament type for an inornate species); A-D denote polymorphic species [A, states 1 + 2; B, states 1 + 3; C, states 2 + 
3; D, states 6 + 7]; BMNH, The Natural History Museum, London]. 
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0 Hormotoma simulatrix (Billings, 1865) 2 5 3 3 

0 Ectomaria adelina (Billings, 1865) 2 5 3 3 

1 Pagodospira cicelia (Billings, 1865) 2 5 3 3 

2 Lophospiraperangulata (Hall, 1847) 2 5 3 3 

3 L. sorrorcula (Billings, 1865) 2 5 3 3 

4 L. rectistriata (Billings, 1865) 2 4 2 3 

5 Pagodospira derwiduii Grabau, 1922 2 5 3 3 

6 Lophospira milleri (Hall in Miller, 1889) 2 4 2 3 

7 Trochonemella knoxvillensis Ulrich in Ulrich and Scofield, 1897 2 5 3 3 

8 Proturritella bicarinata (Koken, 1889) 2 4 2 3 

9 Pagodospira dorothea Grabau, 1922 2 5 3 3 

10 Trochonemella trochonemoides (Ulrich in Ulrich and Scofield, 1897) 2 5 2 3 

11 Lophospira centralis Ulrich and Scofield, 1897 2 3 2 2 

12 Trochonema bellula Ulrich and Scofield, 1897 2 4 2 3 

13 T. eccentrica Ulrich and Scofield, 1897 2 3 2 2 

14 Trochonemella montrealensis Okulitch, 1935 2 4 1 3 

15 Lophospira helicteres (Salter, 1859) 2 3 2 2 

16 Trochonemella n. sp. 2 4 13 

17 Lophospira serrulata (Salter, 1859) 2 5 3 3 

18 Eunema strigillata Salter, 1859 2 3 2 2 

19 Trochonema umbilicata (Hall, 1847) 2 3 2 2 

20 T. canadensis Wilson, 1951 2 3 1 2 

21 T. wilsoni Steele and Sinclair, 1971 2 3 2 2 

22 Lophospira ventricosa (Salter, 1859) 2 5 3 3 

23 Trochonemella notablis (Ulrich in Ulrich and Scofield, 1897) 2 4 1 3 

24 Gyronema pulchella Ulrich and Scofield, 1897 2 3 1 2 

25 G. semicarinata (Ulrich and Scofield, 1897) 2 3 1 2 

26 G. liljevalli Rohr, 1980 2 3 2 2 

27 Trochonema madisonense Ulrich and Scofield, 1897 2 4 2 3 

28 Lophospira burginensis Ulrich and Scofield, 1897 2 5 3 3 

29 L. oweni Ulrich and Scofield, 1897 2 5 3 3 

30 L. concinula Ulrich and Scofield, 1897 2 3 1 2 

31 L. spironema Ulrich and Scofield, 1897 2 4 2 3 

32 Donaldiella decursa (Ulrich and Scofield, 1897) 2 5 3 3 

33 D. conoidea (Ulrich and Scofield, 1897) 2 5 3 3 

34 D. producta (Ulrich and Scofield, 1897) 2 5 3 3 

35 D. Curdsville sp. 2 5 3 3 

36 Lophospira sumnerensis (Safford, 1869) 2 5 3 3 

37 Trochonema salteri Ulrich and Scofield, 1897 2 3 2 2 

38 Ruedemannia humilis (Ulrich and Scofield, 1897) 2 4 2 3 

39 R. lirata (Ulrich and Scofield, 1897) 2 4 2 3 

40 Lophospira tenuistriata Ulrich and Scofield, 1897 2 5 3 3 
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APPENDIX 2. (continued) 


41 L tropidophora Meek, 1873 

2 

5 

3 

3 

1 

1 

2 

1 

2 

2 

2 

2 

2 

2 

1 

42 Donaldiella filosa (Donald, 1902) 

2 

5 

3 

3 

1 

1 

2 

1 

2 

2 

2 

3 

2 

2 

2 

43 D. bowdeni (Safford, 1869) 

2 

5 

3 

3 

1 

1 

2 

1 

3 

2 

2 

3 

2 

2 

2 

44 Lophospira bellicarinata Donald, 1906 

2 

4 

1 

3 

1 

1 

3 

1 

2 

2 

2 

3 

2 

2 

4 

45 L. quadrisulcata Ulrich and Scofield, 1897 

2 

3 

1 

2 

1 

1 

3 

1 

2 

2 

2 

3 

2 

2 

3 

46 “ Schizolopha ” moorei Ulrich and Scofield, 1897 

2 

5 

3 

3 

1 

1 

4 

1 

2 

2 

2 

4 

2 

2 

1 

47 Lophospira aff. serrulate [in Rohr, 1988] 

2 

4 

1 

3 

1 

1 

4 

1 

2 

2 

2 

3 

2 

2 

3 

48 Trochonemella BMNH 36364 

2 

4 

1 

3 

1 

1 

4 

1 

2 

2 

2 

3 

2 

2 

3 

49 T. churkiniRohi, 1988 

2 

4 

2 

3 

1 

2 

4 

1 

2 

2 

2 

3 

1 

. 

- 

50 Trochonemella reusingi Rohr and Blodgett, 1985 

2 

4 

2 

3 

1 

2 

4 

1 

2 

2 

2 

3 

1 

- 

- 

51 Trochonema aff. umbilicata [in Rohr, 1988] 

2 

2 

2 

2 

1 

2 

3 

1 

1 

1 

- 

- 

2 

2 

3 

52 Arjamannia thraivensis (Longstaff, 1924) 

2 

4 

1 

3 

1 

1 

3 

1 

2 

2 

2 

3 

2 

2 

4 

53 Globonema bicarinata (Wenz, 1938) 

2 

1 

1 

2 

1 

1 

3 

1 

1 

1 

- 

- 

2 

2 

4 

54 Lophospira ?borkholmensis (Koken, 1925) 

2 

4 

2 

3 

1 

1 

2 

1 

4 

2 

2 

2 

2 

1 

2 

55 Trochonema pandori Koken, 1925 

2 

4 

2 

3 

1 

2 

3 

1 

1 

1 

- 

- 

2 

2 

3 

56 T. aff. pandori Koken, 1925 

2 

4 

2 

3 

1 

2 

3 

1 

1 

1 

- 

- 

2 

2 

3 

57 Lophospira sedgewicki Donald, 1905 

2 

4 

1 

3 

1 

1 

3 

1 

2 

2 

2 

3 

2 

2 

4 

58 L. Ihyaecinthinsis Foerste, 1924 

2 

4 

1 

3 

1 

1 

3 

1 

2 

2 

2 

3 

2 

2 

4 

59 Longstaffia centervillensis (Foerste, 1923) 

2 

4 

1 

3 

1 

1 

4 

1 

2 

2 

2 

2 

2 

2 

4 

60 Arjamannia cancellata (M’Coy in Sedgwick and M’Coy, 1852) 

2 

4 

1 

3 

1 

1 

3 

1 

2 

2 

2 

2 

2 

2 

4 

61 A. woodlandi (Longstaff, 1924) 

2 

4 

1 

3 

1 

1 

3 

1 

2 

2 

2 

2 

2 

2 

4 

62 Donaldiella trilineata (Foerste, 1923) 

2 

5 

3 

3 

1 

1 

2 

1 

3 

2 

2 

3 

2 

2 

1 

63 Lophospira gotlandica Ulrich in Ulrich and Scofield, 1897 

2 

4 

1 

3 

1 

1 

3 

1 

2 

2 

2 

3 

2 

2 

4 

64 Longstaffia laquetta (Lindstrom, 1884) 

2 

4 

1 

3 

2 

1 

4 

2 

2 

2 

2 

2 

2 

1 

4 

65 Arjamannia inexpectans (Hall and Whitfield, 1872) 

2 

4 

1 

3 

1 

1 

3 

1 

2 

2 

2 

2 

2 

2 

4 

66 Lophospira holmi (Lindstrom, 1884) 

2 

4 

1 

3 

1 

1 

3 

1 

2 

2 

2 

2 

2 

2 

4 

67 Kiviasukkaan nelsonae Peel, 1975b 

2 

4 

1 

3 

1 

1 

3 

1 

2 

2 

2 

3 

2 

2 

3 

68 Lophospira munda (Lindstrom, 1884) 

2 

4 

1 

3 

2 

1 

3 

3 

2 

2 

2 

3 

2 

2 

3 

69 L. imbricate (Lindstrom, 1884) 

2 

4 

1 

3 

2 

1 

4 

4 

2 

2 

2 

3 

2 

2 

3 

70 Trochonema turrita (Lindstrom, 1884) 

1 

- 

- 

- 

2 

1 

4 

1 

1 

1 

- 

- 

2 

2 

3 

71 Eunema kayesi Rohr, 1981 

1 

- 

- 

- 

2 

1 

4 

1 

1 

1 

- 

- 

2 

2 

3 

72 Ruedemannia laevissima (Lindstrom, 1884) 

2 

4 

2 

3 

1 

1 

2 

1 

4 

2 

2 

2 

2 

1 

1 

73 Arjamannia aulangonensis (Peel, 1975a) 

2 

4 

1 

3 

1 

1 

3 

1 

3 

2 

2 

2 

2 

2 

4 

74 Ruedemannia robusta (Lindstrom, 1884) 

2 

4 

2 

3 

1 

1 

2 

1 

4 

2 

2 

2 

2 

1 

1 

75 Longstaffia tubulosa (Lindstrom, 1884) 

2 

4 

1 

3 

2 

1 

4 

3 

2 

2 

2 

2 

2 

1 

4 

76 Longstaffia cyclonema ( Salter, 1873) 

2 

4 

1 

3 

1 

1 

4 

1 

2 

2 

2 

2 

2 

1 

4 

77 Trochonema fatua (Whiteaves, 1895) 

2 

1 

2 

2 

1 

1 

3 

1 

1 

1 

- 

- 

2 

2 

3 

78 Loxoplocus soluta (Whiteaves, 1884) 

2 

4 

1 

3 

1 

1 

3 

1 

2 

2 

2 

3 

2 

2 

2 

79 Lophospira cochleata (Lindstrom, 1884) 

2 

4 

1 

3 

2 

1 

4 

3 

2 

2 

2 

3 

2 

2 

3 

80 Eunema muricata (Lindstrom, 1884) 

1 

- 

- 

- 

2 

1 

4 

1 

1 

1 

- 

- 

2 

2 

3 

81 Ruedemannia subrobusta (Pemer, 1907) 

2 

4 

2 

3 

1 

1 

2 

1 

4 

2 

2 

2 

2 

1 

1 

82 Ptychozone aberrans Pemer, 1907 

2 

4 

2 

3 

1 

1 

2 
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1 
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5 
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5 

5 

5 

5 
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6 

7 

8 

9 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

o 
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0 H. simulatrix 

? 
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2 
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- 

- 
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1 

0 E. adelina 

? 

2 

1 

- 

- 

1 

1 

- 

- 

1 

l 

1 

- 

3 

2 

5 

2 

3 

3 

1 

5 

3 

1 

- 
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- 

- 

- 

- 

? 

? 
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1 

3 

2 

2 
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1 

1 
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1 P. cicelia 

1 

2 

1 

- 

- 

1 

2 

- 

- 

- 

1 

2 

1 

3 

3 

5 

2 

3 

3 

1 

5 

3 

1 

- 

1 

- 

- 

- 

- 

? 

7 

? 

1 

3 

2 

2 

A 

1 

1 

1 

2 L. perangulata 

1 

2 

1 

- 

- 

1 

2 

- 

- 

- 

l 

2 

1 

3 

3 

5 

2 

3 
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3 

1 
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1 

- 

- 

- 

- 

? 

? 
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1 
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2 
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A 

1 

1 

1 

3 L. sorrorcula 

l 

2 

1 

- 
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1 

2 

- 

- 

- 

l 

2 

1 

3 

3 

4 

2 

3 

3 

1 

5 

3 

1 
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1 

- 

- 

- 

- 

7 

? 

? 

1 

3 

2 

2 

1 

1 

1 

1 

4 L. rectistriata 

1 

2 

1 

- 

- 

1 

2 

- 

- 

- 

1 

2 

1 

3 

3 

4 

2 

3 

3 

1 

5 

3 

1 

- 

1 

- 

- 

- 

- 

? 
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7 

1 

3 

2 
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1 

1 

5 P. denviduii 

1 

3 

1 

- 

- 

1 

2 

- 

- 

- 

1 

2 

1 

3 

3 

5 

2 

3 

3 

1 

5 

3 

1 

- 

1 

- 

- 

- 

- 

? 

7 

7 

1 

3 

2 

1 

- 

- 

1 

1 

6 L. milleri 

1 

2 

1 

- 

- 

1 

2 

- 

- 

- 

1 

2 

1 

3 

3 

4 

3 

3 

4 

1 

5 

3 

1 

- 

2 

2 

3 

1 

1 

1 

1 

1 

1 

3 

1 

2 

2 

2 

1 

1 

1 T. knoxvillensis 

1 

2 

1 

- 
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2 

- 
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2 

2 

6 
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A 
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1 

8 P. bicarinata 

1 

2 

1 

- 

- 
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1 

2 
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3 

3 

4 

3 

3 

4 

1 
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9 P. dorothea 
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3 
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3 

1 
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1 
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1 

3 

2 
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- 

- 

1 
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10 T. trochonemoides 

1 

2 

1 

- 

- 

1 

2 

- 
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- 

1 

2 

2 
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2 
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2 
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2 
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1 

1 

1 
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1 

3 
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2 

A 

1 

1 

1 

11 L. centralis 

1 

2 

1 

- 

- 

4 

2 

- 

- 

- 

1 

2 

1 

3 

2 

3 

3 

3 

4 

1 

5 

3 

1 

- 

2 

2 

3 

1 

1 

1 

1 

1 

1 

3 

1 

2 

2 

2 

1 

1 

12 T. be Hula 

1 

2 

1 

- 

- 

1 

2 

- 

- 

- 

1 

2 

2 

7 

2 

3 

2 

7 

7 

1 

7 

5 

2 

2 

2 

2 

3 

1 

2 

2 

l 

1 

1 

3 

1 

2 

2 

2 

1 

1 

13 T. eccentrica 

1 

2 

1 

- 

- 
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- 

- 
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3 
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1 

14 T. montrealensis 

l 
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- 
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- 
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7 
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3 
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1 

1 

15 L. he lie teres 

1 

2 

1 

- 

- 

4 

2 

- 

- 

- 

l 

2 

1 

3 

2 

3 

3 

3 

4 

1 
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2 
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2 

3 

3 
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2 

1 

1 

1 

1 

3 

1 

2 

2 

2 

1 

1 

16 T. n. sp. 

l 

2 

2 
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2 
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2 

2 

1 

1 

3 
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1 

l 

17 L. serrulata 

2 

2 

1 

- 

- 

1 

? 

2 

2 

- 

1 

2 

1 

3 

3 

5 

2 

3 

3 

1 

5 

3 

2 

1 

2 

3 

3 

1 

1 

1 

1 

1 

1 

3 

1 

2 

2 

2 

1 

1 

18 E. strigillata 

1 

2 

1 

- 

- 
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2 

- 

- 

- 
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2 

2 
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2 

3 

3 

3 

4 

1 

5 

3 

2 

1 

2 

3 

3 

1 

2 

1 

1 

1 

1 

3 

1 

2 

2 

2 

1 

1 

19 T. umbilicata 

1 
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- 
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- 

- 
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2 
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7 
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2 

7 

7 

1 

7 
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2 

1 

1 

l 
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2 

2 

1 

l 

20 T. canadensis 

1 

2 
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- 

- 
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- 
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2 

2 
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2 

7 

7 

1 

7 
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1 

3 
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21 T. wilsoni 
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- 

- 
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7 

2 
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2 

7 

7 

1 

7 

6 
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2 
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2 
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1 

1 

3 

1 

2 

2 

2 

1 

l 

22 L. ventricosa 

1 
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3 
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3 

1 
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l 

23 T. notablis 
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1 

24 G. pulchella 
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- 
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2 

3 

3 

3 

4 

1 
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25 G. semicarinata 
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26 G. liljevalli 
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- 

1 

- 

. 

- 

. 

1 

3 

L. sorrorcula 

1 

2 

1 

- 

- 

- 

? 

3 

3 

2 

1 

1 

1 

3 

2 

1 

- 

1 

- 

1 

- 

1 

- 

3 

6 

6 

3 

2 

2 

- 

1 

1 

- 

- 

1 

- 

_ 

_ 

_ 

2 

4 

L. rectistriata 

1 

2 

1 

- 

- 

- 

? 

3 

3 

2 

1 

1 

1 

3 

2 

1 

- 

1 

- 

1 

- 

1 

_ 

3 

6 

6 

3 

2 

2 



1 



1 





2 

5 

P. derwiduii 

- 

? 

1 

- 

- 

- 

? 

2 

3 

2 

1 

1 

1 

3 

? 

1 

- 

1 

- 

1 

_ 

1 

- 

2 

5 

8 

3 

2 

C 

1 

1 

1 

_ 

_ 

1 

_ 

_ 

_ 

_ 

1 

6 

L. milleri 

1 

2 

1 

- 

- 

- 

? 

4 

3 

2 

1 

2 

1 

3 

2 

1 

- 

1 

- 

1 

- 

1 

- 

3 

6 

6 

3 

2 

2 



1 



1 





1 

7 

T. knoxvillensis 

1 

2 

2 

1 

2 

1 

2 

4 

2 

2 

1 

4 

1 

4 

1 

2 

3 

1 

. 

1 

_ 

2 

2 

4 

4 

5 

3 

1 

3 

1 

1 

1 

_ 

_ 

1 

_ 


_ 

_ 

2 

8 

P. bicarinata 

1 

2 

1 

- 

- 

- 

? 

4 

3 

2 

1 

2 

1 

3 

2 

1 

- 

1 

- 

1 

- 

1 

- 

3 

6 

6 

3 

2 

2 

- 

- 

2 

3 

3 

2 

4 

3 

1 

1 

1 


(continued) 








? 

1 

2 

1 

1 

2 

2 

2 

2 

2 

1 

1 

1 

2 

2 

2 

2 

1 

1 

3 

4 

2 

2 

2 

2 

2 

2 

4 

1 

2 

2 

3 

4 

2 

3 

2 

2 

9 

2 

2 

2 

2 

1 

2 

1 

? 

1 

1 

2 

2 

2 

2 

2 

2 

3 

2 

2 

2 

3 

2 

2 

1 

1 

1 

? 

1 

2 

2 

1 

? 

3 

1 

1 
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APPENDIX 3 
STRATIGRAPHIC DATA 

Range and sampling data for lophospiroid species. Separate horizon scales are used for each realm when calculating the likelihood of gaps (Wagner, 1995a). 
Stratigraphic debt and discrete ranges were calculated using the scales on Fig. 1, where Cassinian = Middle Arenig, Fennian = Late Arenig, Llanvim = Early 
Llanvim, Llandeilo = Late Llanvim, Ashbian-Black Riveran = Early Caradoc, Rocklandian-Kirkfieldian = Middle 1 Caradoc, Shermanian-Edenian = 
Middle2 Caradoc, Maysvillian-Richmondian = Ashgill, Rhuddanian-Aeronian = Early Llandovery, Telychian = Late Llandovery, Sheinwoodian = Early 
Wenlock, Homerian = Late Wenlock, Gorstian = Early Ludlow, and Ludfordian = Late Ludlow. [FKA, oldest sampled horizon, counted as number of older 
lophospiroid horizons known up to that point (see Wagner, 1995a); FSA, first stage from which a species is known; H, number of horizons from which a 
species is known; LKA, latest sampled horizons; LSA, last stage from which a species is known; Realm, biogeographic province: Baltica (Balt), 
Gondwana (Gond), Laurentia (Laur), Toquima-Tablehead (ToqTab)]. 


No. Species 

H 

FKA 

LKA 

Realm 

FSA 

LSA 

1 

Pagodospira cicelia 

14 

1 

57 

Laur 

Cassinian 

Llanvim 

2 

Lophospira perangulata 

85 

1 

256 

Laur 

Cassinian 

Richmondian 

3 

L. sorrorcula 

16 

18 

57 

ToqTab 

Fennian 

Llandeilo 

4 

L. rectistriata 

6 

18 

44 

Laur 

Fennian 

Llandeilo 

5 

Pagodospira derwiduii 

7 

18 

27 

ToqTab 

Fennian 

Llanvim 

6 

Lophospira milleri 

49 

18 

256 

Laur 

Llanvim 

Richmondian 

7 

Trochonemella knoxvillensis 

7 

18 

44 

Laur 

Llanvim 

Llandeilan 

8 

Proxurritella bicarinata 

19 

22 

48 

Balt 

Llandeilan 

Llandeilan 

9 

Pagodospira dorothea 

1 

22 

44 

ToqTab 

Llandeilo 

Llandeilan 

10 

Trochonemella trochonemoides 

5 

22 

46 

Laur 

Llandeilo 

Ashby 

11 

Lophospira centralis 

35 

45 

220 

Laur 

Ashby 

Rocklandian 

12 

Trochonema bellula 

1 

45 

46 

Laur 

Ashby 

Black Riveran 

13 

T. eccentrica 

1 

45 

47 

Laur 

BlackRiveran 

Black Riveran 

14 

Trochonemella montrealensis 

8 

47 

91 

Laur 

BlackRiveran 

Black Riveran 

15 

Lophospira helicteres 

16 

47 

131 

Laur 

BlackRiveran 

Rocklandian 

16 

Trochonemella n. sp. 

3 

47 

91 

Laur 

BlackRiveran 

Black Riveran 

17 

Lophospira serrulata 

16 

47 

137 

Laur 

BlackRiveran 

Kirkfieldian 

18 

Eunema strigillata 

5 

47 

131 

Laur 

BlackRiveran 

Black Riveran 

19 

Trochonema umbilicata 

31 

47 

203 

Laur 

BlackRiveran 

Shermanian 

20 

T. canadensis 

3 

47 

91 

Laur 

BlackRiveran 

Shermanian 

21 

T. wilsoni 

1 

47 

58 

Laur 

BlackRiveran 

Shermanian 

22 

L. ventricosa 

13 

59 

208 

Laur 

BlackRiveran 

Kirkfieldian 

23 

Trochonemella notablis 

4 

92 

131 

Laur 

BlackRiveran 

Black Riveran 

24 

Gyronema semicarinata 

6 

59 

131 

Laur 

Rocklandian 

Rocklandian 

25 

G. pulchella 

4 

92 

137 

Laur 

BlackRiveran 

Rocklandian 

26 

G. liljevalli 

1 

97 

107 

ToqTab 

Rocklandian 

Rocklandian 

27 

Trochonema madisonense 

7 

97 

256 

Laur 

Kirkfieldian 

Richmondian 

28 

Lophospira burginensis 

13 

97 

209 

Laur 

Rocklandian 

Edenian 

29 

L. oweni 

11 

119 

209 

Laur 

Rocklandian 

Edenian 

30 

L. cone inula 

4 

132 

165 

Laur 

Rocklandian 

Rocklandian 

31 

L. spironema 

2 

132 

137 

Laur 

Rocklandian 

Rocklandian 

32 

Donaldiella decursa 

1 

155 

184 

Laur 

Kirkfieldian 

Shermanian 

33 

D. conoidea 

1 

155 

165 

Laur 

Kirkfieldian 

Shermanian 

34 

D. producta 

1 

155 

165 

Laur 

Kirkfieldian 

Shermanian 

35 

D. Curdsville sp. 

1 

187 

208 

Laur 

Kirkfieldian 

Shermanian 

36 

Lophospira sumnerensis 

15 

138 

208 

Laur 

Shermanian 

Richmondian 

37 

Trochonema salteri 

8 

138 

265 

Laur 

Shermanian 

Richmondian 

38 

Ruedemannia humilis 

3 

155 

184 

Laur 

Edenian 

Richmondian 

39 

R. lirata 

6 

187 

253 

Laur 

Edenian 

Richmondian 

40 

Lophospira tenuistriata 

3 

187 

220 

Laur 

Edenian 

Richmondian 

41 

L. tropidophora 

29 

187 

256 

Laur 

Edenian 

Richmondian 

42 

D. filosa 

1 

187 

208 

Laur 

Edenian 

Edenian 

43 

D. bowdeni 

21 

187 

256 

Laur 

Edenian 

Richmondian 

44 

Lophospira bellicarinata 

2 

210 

256 

Laur 

Maysvillian 

Richmondian 

45 

L. quadrisulcata 

6 

210 

220 

Laur 

Maysvillian 

Maysvillian 

46 

L. moorei 

7 

221 

256 

Laur 

Richmondian 

Richmondian 

47 

L. aff. serrulata 

3 

221 

256 

Laur 

Richmondian 

Richmondian 

48 

Trochonemella BMNH 36364 

1 

221 

256 

Laur 

Richmondian 

Richmondian 

49 

T. churkini 

3 

221 

256 

Laur 

Richmondian 

Richmondian 


(continued) 
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APPENDIX 3. (continued) 

50 T. re us in gi 

3 

221 

256 

Laur 

Richmondian 

Richmondian 

51 

Trochonema aff. umbilicata 

3 

221 

256 

Laur 

Richmondian 

Richmondian 

52 

Arjamannia sybellina 

5 

221 

256 

Laur 

Richmondian 

Richmondian 

53 

Globonema bicarinata 

2 

223 

239 

Balt 

Richmondian 

Richmondian 

54 

Lophospira Iborkholmensis 

2 

223 

239 

Balt 

Richmondian 

Richmondian 

55 

Trochonema pandori 

2 

223 

239 

Balt 

Richmondian 

Richmondian 

56 

T. aff. pandori 

1 

223 

239 

Balt 

Richmondian 

Richmondian 

57 

Lophospira sedgewicki 

5 

223 

273 

Laur 

Richmondian 

Telychian 

58 

L. 'Ihyacinthensis 

2 

223 

256 

Laur 

Richmondian 

Richmondian 

59 

Longstaffia centervillensis 

3 

257 

273 

Laur 

Rhuddanian 

Telychian 

60 

Arjamannia cancellata 

1 

257 

261 

Laur 

Rhuddanian 

Telychian 

61 

A. woodlandi 

5 

257 

273 

Laur 

Rhuddanian 

Telychian 

62 

D. trilineata 

3 

257 

265 

Laur 

Aeronian 

Aeronian 

63 

Lophospira gotlandica 

30 

257 

337 

Laur 

Aeronian 

Ludfordian 

64 

Longstaffia laquetta 

6 

257 

301 

Laur 

Rhuddanian 

Homerian 

65 

A. inexpectans 

3 

262 

273 

Laur 

Telychian 

Telychian 

66 

Lophospira holmi 

1 

275 

277 

Laur 

Sheinwoodian 

Sheinwoodian 

67 

Kiviasukkaan nelsonae 

1 

275 

277 

Laur 

Sheinwoodian 

Sheinwoodian 

68 

Lophospira munda 

1 

282 

294 

Laur 

Sheinwoodian 

Homerian 

69 

L. imbricata 

7 

278 

330 

Laur 

Sheinwoodian 

Gorstian 

70 

Eunema turrita 

5 

278 

308 

Laur 

Sheinwoodian 

Homerian 

71 

E. kayesi 

1 

282 

294 

Laur 

Sheinwoodian 

Gorstian 

72 

Ruedemannia laevissima 

1 

295 

301 

Laur 

Homerian 

Homerian 

73 

Arjamannia aulongensis 

1 

295 

301 

Laur 

Homerian 

Whitwellian 

74 

Ruedemannia robusta 

7 

282 

330 

Laur 

Homerian 

Gorstian 

75 

Longstaffia tubulosa 

3 

295 

330 

Laur 

Homerian 

Gorstian 

76 

L. cyclonema 

5 

302 

308 

Laur 

Homerian 

Homerian 

77 

Trochonema fatua 

3 

309 

330 

Laur 

Gorstian 

Gorstian 

78 

Loxoplocus soluta 

4 

309 

330 

Laur 

Gorstian 

Gorstian 

79 

Lophospira cochleata 

3 

309 

330 

Laur 

Gorstian 

Gorstian 

80 

Trochonema muricata 

1 

331 

332 

Laur 

Ludfordian 

Ludfordian 

81 

Ruedemannia subrobusta 

3 

331 

337 

Gond 

Ludfordian 

Ludfordian 

82 

Ptychozone aberrans 

2 

331 

337 

Gond 

Ludfordian 

Pridolf 


APPENDIX 4 


SYSTEMATIC PALEONTOLOGY 

In the interest of conserving space, taxonomic revisions are pre¬ 
sented only for new genera and for those requiring rediagnoses. Genera 
merely redefined but retaining traditional diagnoses are reclassified as per 
Fig. 14. 

SUPERFAMILY LOPHOSPIROIDEA NOM. TRANS. WENZ, 1938 
FAMILY LOPHOSPIRIDAE WENZ, 1938 

Diagnosis: Lophospiroids primitively featuring a trilineate sinus 
keel, a strong left ramp carina, and a deep sinus that curves back to the 
sinus keel. All of these features change within the family in at least one 
subclade. 

Included Genera: Lophospira, Donaldiella, Loxoplocus, 
Eunema, Gyronema, Ruedemannia, Arjamannia, Paupospira gen. nov., 
Frodospira gen. nov. 

Genus Eunema Salter, 1859 

Diagnosis: Prominent, sharp right and left ramp carina. Strong 
medial keel usually bordered by peripheral lira (although not in types 
species). Shallow sinus trending nearly straight into sinus keel. 
Increasing translation during final whorls, sometimes resulting in disjunct 
coiling. 

Type Species: Eunema strigillata Salter, 1859. 


Additional Species: "Lophospira ” centralis, “L. ” helicteres, 
“L. ” concinula Ulrich and Scofield, 1897, “L. ” quadrisulcata. 

Discussion. Knight et al. (1960) considered Eunema to be a poly- 
phyletic subgenus of Trochonema, meant to describe high-spired variants 
of that genus. However, phylogenetic analyses indicate that the synapo- 
morphies of Eunema and Trochonema are parallelisms. Knight et al. 
(1960) acknowledged the polyphyletic nature of the Trochonematidae by 
noting that the diagnostic characters ( e. g. shallow sinus and a medial keel) 
appeared to be convergent among many forms. This analysis takes their 
conclusions a step further by positing that the form is polyphyletic among 
lophospiroids. The much more restrictive definition presented here 
excludes several species previously assigned to Eunema , which have been 
reassigned to genera such as Globonema. 

First-known appearance (FKA): “ Lophospira ” centralis: 
Murfreesboro Limestone (Early Caradoc [Ashbyan]). 

Last-known appearance (LKA): “Lophospira ” quadrisulcata: 
Maquoketa Formation (Early Ashgill [Maysvillian]). 

Genus Loxoplocus Fischer, 1885 

Syn. Kiviasukkaan Peel, 1975b 

Diagnosis: Sharp and strong medial keel that hooks slightly 
adapically, typically bordered by sharp but much weaker peripheral lira 
(although seemingly not on the type species). Narrow, shallow sinus, 
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curving back towards sinus keel. Dull but strong left ramp carina, with 
right ramp concave and featuring a sharp carina at the suture, which cre¬ 
ates an acute sutural margin. Increasing translation during final whorls, 
sometimes resulting in disjunct coiling (with the type species showing dis¬ 
junct coiling throughout life). 

Type Species: “Murchisonia " soluta Salter, 1859. 

Additional Species: “Lophospira ” sedgewicki Donald, 1906, 
“L. ” Ibellicarinata Donald, 1906, “L.” gotlandica Ulrich and Scofield, 
1897, “Kiviasukkaan" nelsonae Peel, 1975b. 

Discussion: Knight et al. (1960) originally defined Loxoplocus 
as a broad genus that included Lophospira and other subgenera, an inter¬ 
pretation rejected by more recent authors (e. g. Tofel and Bretsky, 1987). 
It now is restricted to a paraclade of lophospiroids whose most prominent 
features include increasing translation over ontogeny, a very acute suture, 
and a very concave right ramp. 

FKA: “ Lophospira” sedgewicki: “Starfish Beds,” Girvan 
District (Late Ashgill [Rawtheyan]). 

LKA: “ Lophospira” gotlandica: Kopanina Formation (Late 
Ludlow [Ludfordian]); 

Genus Paupospira gen. nov. 

Fig. 14 

Diagnosis: Thick columella, filling the entire umbilicus and pro¬ 
ducing a shovel-like siphon in extreme cases. Deep and wide sinus curv¬ 
ing back strongly to sinus keel. Sinus keel trilineate on early whorls, but 
with trilineations often vague and weak on adult whorls. 

Type Species: Lophospira oweni Ulrich and Scofield, 1897. 

Additional Species: “Lophospira ” burginensis , “L. ” tenuistriata 
Ulrich and Scofield, 1897, “L. ” sumnerensis Ulrich and Scofield, 1897, 
“L. ” tropidophora Meek, 1873, “ Schizolopha ” moorei, “Donaldiella ” 
bowdeni (Safford, 1869), “ Hormotoma ” trilineata Foerste, 1923. 

Etymology: For David Swofford’s computer program PAUP 
(“Phylogenetic Analysis Using Parsimony”; Swofford, 1998), which 
helped to diagnose and define this clade. 

Discussion: Paupospira matches a monophyletic group original¬ 
ly hypothesized by Ulrich and Scofield (1897). A new genus is erected to 
recognize that group. The most distinguishing feature is the extremely 
thick columella, which often preserves as a core without the rest of the 
shell. Cladistic depictions of relationship highlighting the key features are 
given in Fig. 14. 

FKA: “ Lophospira ” burginensis: Leray beds, Rockland 


Formation (Middle Caradoc [Rocklandianj). 

LKA: “ Hormotoma ” trilineata: Saugh Hill Group (Middle 
Llandovery [Aeronian]). 

Genus Frodospira gen. nov. 

Diagnosis: Very small shells, with strong ornament that leaves 
imbrications on the growth lines, a strongly adapically curved medial keel, 
and highly tangential apertures. 

Type Species: “Murchisonia " imbricata Lindstrom, 1884. 

Additional Species: “ Lophospira" munda (Lindstrom, 1884), 
“L. ” cochleata (Lindstrom, 1884). 

Etymology: After J. R. R. Tolkein’s character from Lord of the 
Rings , reflecting the unusually small size of the known species. 
Hypotheses about extraneous epidpodial tentacles are purely speculative. 

Discussion: This closely matches Ulrich and Scofield (1897: 
963) imbricata subsection of Lophospira. The genus is remarkable in that 
all known species are known only from very small shells. As these shells 
are found in the same assemblages as large gastropods, it appears to repre¬ 
sent a character of the genus rather than a taphonomic artifact. The illus¬ 
trations of these species, provided by Lindstrom (1884) are very accurate 
and convey the general characters very well. 

FKA: “ Lophospira ” munda: C beds, Visby Formation (Early 
Wenlock [Sheinwoodian]). 

LKA: “Lophospira” imbricata: Upper Hemse Beds (Early 
Ludlow [Gorstian]). 

FAMILY TROCHONEMATIDAE ULRICH 
AND SCOFIELD, 1897 

Diagnosis: Trochiform, widely umbilicate lophospiroids. Early 
forms retain wide, and deep sinus curving back to a trilineate selenizone, 
but later forms possess either the medial lira only or the peripheral lira 
only. Sharp left and right ramp carina with strong umbilical and sutural 
carina found on most species. 

Included Genera: Trochonemella, Trochonema, Globonema. 

Discussion: Knight et al.’s (1960) reduction of the 
Trochonematidae is advanced here. Preliminary analyses (Wagner, 
unpub. data) indicate that Devonian genera assigned to the 
Trochonematidae such as Trochonemopsis belong to the Euomphalinae. 
At this time, it is not clear if any post-Silurian genera assigned to this fam¬ 
ily belong here. 


